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AUTHOR'S  SUMMARY  '  '  ^ 

Quartz  oscillators,  used  for  aerospace  applications,  are  subjected  to  very 
severe  environmental  conditions,  particularly  as  far  as  accelerations  and  vibra¬ 
tions,  which  can  be  periodic  or  random,  are  concerned.  The  main  effect  of  these 
disturbances  is  to  modify  the  natural  frequency  of  the  resonator  and  thus  to 
cause  frequency  instabilities  in  the  oscillator. 


In  this  paper  the  quartz  resonator  is  studied  both  by  itself,  as  a  passive 
element  placed  in  a  transmitting  circuit,  and  installed  in  an  oscillator,  which 
is  its  normal  condition  of  usage. 


The  accelerations  to  which  the  resonators  were  subjected  were  steady  state, 
sinusoidal  or  random.  The  orientation  of  the  acceleration  relative  to  the  quartz 
crystallographic  axis,  as  well  as  its  direction,  are  of  prime  importance. 


The  same  tests  were  performed  on  oscillators  containing  the  resonators 
previously  studied.  The  results  confirm  the  part  played  by  the  resonator  in  the 
degradation  of  the  stability  of  a  quartz  oscillator  subjected  to  external 
disturbances . 
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INTRODUCTION 


The  primary  standard  used  at  the  present  time  for  frequency  measurement  is 
the  caesium  clock.  In  general,  atomic  oscillators  are  used  as  primary  frequency 
standards.  However,  recent  progress  into  the  long-term  stability  of  quartz 
oscillators  has  given  them  performances  which  approach,  closely,  those  of  quantum 
oscillators. 

Time-frequency  systems  have  numerous  applications  such  as  for  aerial  and 
maritime  navigation,  anticollision  systems,  flights  in  formation,  navigation 
satellites,  all  of  which  require  frequency  standards  to  be  fitted  to  terrestrial 
or  spacecraft.  The  necessary  conditions  of  volume,  weight  and,  also,  price, 
allied  to  their  excellent  characteristics,  lead  to  a  preference  for  quartz 
oscillators  rather  than  atomic  oscillators. 

Aerospace  applications  of  quartz  oscillators  necessitate  the  study  of  the 
influence  of  the  environment  on  oscillator  stability.  The  main  disturbances, 
accelerations  and  vibrations,  mod’fy  the  oscillator  characteristics  by  their 
effects  on  the  quartz  resonator. 

This  paper  begins  by  referring  to  the  general,  nonlinear  theory  of  an  elastic 
body  subjected  to  a  predeformation.  This  theory  allows  us  to  give  an  interpreta¬ 
tion  of  the  experimental  phenomena  observed. 

The  first  part  of  the  experimental  study  deals  with  the  influence  of 
steady-state,  sinusoidal  or  random  accelerations  applied  to  the  resonator  alone. 

In  all  cases  a  disturbance  of  the  resonator  frequency  was  seen.  The  problem  of 
crystal  suspension  is  approached  from  the  very  particular  case  of  two  supports. 

The  second  part  is  devoted  to  the  use  of  the  resonator  in  an  oscillation 
loop.  The  oscillator  so  formed  is  subjected  to  the  same  disturbances  as  was  the 
resonator . 

All  the  results  confirm  the  role  of  the  resonator  in  the  degradation  of  the 
stability  of  a  quartz  oscillator  subjected  to  external  disturbances. 


Part  I  -  STUDY  OF  THE  RESONATOR  SUBJECTED  TO  ENVIRONMENTAL  CONDITIONS 


Chapter  1  Theoretical  aspects  of  the  phenomena 


The  acceleration  applied  to  the  resonator  caused  prestreses  and  predeforma¬ 
tions  due  to  the  internal  forces  and  surface  stresses  (these  latter  being  due  to 
the  fixing  points).  In  the  case  of  a  steady  acceleration  these  prestresses  are 
purely  static.  From  the  mathematical  point  of  view  this  is  almost  true,  also. 
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the  sinusoidal  or  random  accelerations  since  their  maximum  frequencies  (5000  Hz) 
are  still  small  with  respect  to  the  natural  frequency  of  the  resonator  (5  MHz). 

This  problem  is  analogous  to  that  of  the  application  of  an  external,  static 
stress  to  the  crystal,  as  mentioned  by  Lee  ^  or  to  the  problem  of  thermal 

stresses. 

In  this  study  we  use  the  nonlinear  relationships  for  quartz  developed  by 
2 

K.N.  Thurston"". 

Three  states  of  the  crystal  will  be  distinguished:- 

the  natural  state  where  the  crystal  is  not  subjected  to  any  force, 
the  initial  state  where  the  crystal  is  subjected  to  a  static  prestress, 
the  final  state  where  a  dynamic  deformation,  assumed  to  be  of  infinitely- 
small  amplitude,  is  superimposed  on  the  static  deformation. 

The  coordinates  of  a  point  in  the  crystal  in  its  three  states  will  be 
defined  by  a.,  x^  and  ,  and  the  corresponding  densities  by  p^,  p  and  p  . 

The  static  displacement  can  be  defined  as:- 


u .  =  X .  -  a .  . 

1  11 


The  dynamic  displacement  as:- 

u.  =  X.  -  X. 
1  11 


and  the  total  displacement  as:- 


U.  =  X.  -  a.  . 
1  11 


In  the  final  state  the  equilibrium  equation  can  be  written  as:- 


P 


Fi«,) 


+ 


3T.  . 

_ Li 

3X. 

J 


where  F.  represents  the  internal  forces 

and  T..  the  total  stresses, 
ij 

In  the  natural  state  this  equation  can  be  written  as:- 


P 


0 


JF.(X^) 


(1) 

(2) 

(3) 


(A) 


(5) 
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whiTo  J  represonl  s  llio  Jacobii'ii  iiet»’niiiiiaiit 


J  = 


'  U 


li.  till'  samo  way  wt‘  can  also  introJiici' 


J  =  ^ 


It’) 


and  I’.j  is  ilu'  t  irst  I’iola-Kirchliol  1  tensor 


.>X. 


'■>,  \if 


(7) 


^ij  ‘‘’PresiMit  s  the  Ihenuodynani  c  tensions  which  can  he  related  to  the  stresses 

I  .  .  bv  t  he  I'ouat  i  on 
tj 


.*a  .la 

t  =  J  — ^  r  .  , 

''x.  '' 


i8) 


Ihe  therr.iodynamic  tensions  have  considerable  magnitudes  since  it  is  through 
them  that  the  noniial  elastic  coel  1  i  c  i  .-nt  s  are  introduced. 


‘  J/g 


f) 


where  t'  is  tlie  initial  energy  j’er  unit  volume 
and 

t  .  •  =  C.  .  n,  +  U: .  .  n,  ,n  + 

ij  ijk>k>  i.)k»  mn  k  i  mn 


(.10) 


where  represents  the  deformat  ions . 

The  first  I’iola-Ki  rchhof  f  tensor  1’..  is  a  function  of  the  gradients 
enly.  It  can  also  be  written  in  the  form  ol  a  series  development 
limited  to  the  first  two  terms:- 
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+ 


’a 

m 


,11' _ 


U  0 


8 


wluTi'  r.  .  is  till'  valui'  ot 
1  J 

It  Wl'  1  f  t 


r.j  i'X|ui'SSt'ii  ill  till'  initial  state. 


the  i'i[u  i  1  i  In  i  uiii  i  iin.i  t  i  an 


t  ''K  ean  be  wi  i  i  I  en 


u  . 
i 


JK.  t.\  ) 
1  1 


JK. tx  ) 

i  I 


and  by  putting  eiiuation  (7)  into  equat ion  tIJ)  we  have:- 


A. 


i  k  i  m 


‘X.  ,ix 

■  .  t  .  +  - -  — C  . 

ik  jm  da  d,i  j  pniq 
1’  4 


U2) 


(1  }) 


U-i) 


i  the  Kriineekei  symbol. 

Ihe  equilibtiuui  equation  t.  1  1)  is  expressed  in  the  natural  state.  It  is 
at  interest  to  express  it  in  the  initial  state  as:- 


,n 


U.  _  / - ,lu 

“T"  ~  ~  1' •  t  M  “  K  .  t  X  )  +  -7^ —  (  i  sk  r  ^ 

it'  I  '  eX  \  dx 


r  / 


where  we  have  let:- 


15. 


i  skr 


I'i..  T 
ik  sr 


i  skr 


and 


^  ilrnu^uv  - 

5'  .  UV 

J  1’>"‘1 


'i  skr 


?x. 

dx 

dx,  dx 

1 

s 

_ k  _ r 

da  . 

da 

da  da 

J 

P 

m  q 

t: . 


U3) 


Ut’) 


U7) 


In  the  ease  under  eonsideration,  we  can  consider  that  the  accelerations 
.ire  independent  of  tlie  deformations  and  hence: 


K.  (,X  ) 

i  { 

f  .  ix  , ) 

i  t 


pl' 


fI8) 

fid)  i;r 

18d7 


V 


where  I' 
The 


is  the  acceleration  applied, 
equilibrium  equation  then  reduces 


to:- 


■) 


A  particularly 
case  the  coefficient 
can  be  written  as:- 


simple  case  is  that  of  homogeneous 


B.  ,  is  independent  of  x  and 
tskr  s 


de format  ion.  In  this 
the  equilibrium  equation 


■) 


3t 


) 

3  u, 


iskr  3x  3x 
s  r 


(21) 


This  equation  is  analagous  to  that  for  a  linear  system  with 
The  coefficients  B.  ,  simplv  replace  the  elastic  coefficients 
advantage  is  that  this  equation  is  of  a  general  form  and  can  thus 
to  any  elastic  model  subjected  to  prestresses  due  to  acceleration 
forces  or  to  temperature  gradients  . 


predef ormat ion . 

C .  ,  .  The 

iskr 

be  applied 
to  external 


It  is  now  possible  to  consider  the  problem  in  two  parts.  The  first  part 
is  the  study  of  the  static  deformation  alone.  Since  the  amplitude  of  the  dynamic 
deformation  is  very  small  (.and  this  is  always  true,  at  least  in  part)  it  can  be 
considered  to  have  no  influence  on  the  static  deformation.  This  latter  can  then 
be  studied  using  a  two-  or  three-dimensional  linear  model  as  needed. 

The  second  part  is  the  study  of  the  dynamic  deformation  itself.  This  can 
be  treated  with  the  aid  of  a  one-dimensional  model  using  equations  (20)  or  (21) 
depending  upon  whether  the  static  deformation  is  homogeneous  or  not. 

Let  us  consider  the  case  of  a  quartz  plate  cut  along  AT  vibrating  in 
shear  across  the  thickness  and  let  us  assume  the  static  deformation  to  be 
homogeneous.  Equation  (21)  reduces  to;- 


1 


(22) 


The  coefficient  B  ^  can  be  expressed  by  using  equations  (16)  and  (17) 


as :  * 


J  ’^2  ^66 


*  2Sj  .  2S2  . 


c  s+c  s+c  s+c  s 

661  I  662  2  663  3  664  4 


'66 


+  2 


^66  ^^3 


(23) 


where  anti  define  the  deformations  and  static  forces  in  single  index 

notation. 

The  vibration  frequency  can  be  written  as:- 


When  there  is  no  prestress,  the  frequency  is  equal  to:- 


f 


0 


I  /^66 

V  '’o 


We  also  have:- 


e 


^2^ 


from  which  we  obtain:- 


f 


J  B 


12  12 


(1  -  2S2) 
66 


(24) 


(25) 


(26) 


(27) 


The  resonators  used  have  electrodes  fitted  at  the  centres  of  the  plates. 

The  plate  can  then  be  considered  to  vibrate  mainly  about  its  centre,  and  of  the 

distribution  of  forces  and  static  deformations,  only  their  values  at  the  centre 

need  be  considered.  It  can  be  shown,  easily,  that  if  these  deformations  are  due 

to  internal  forces  and  if  the  resonator  is  symmetrical  with  respect  to  its 

principle  planes  (Fig  6),  the  deformations  S,,  S  and  S 
-  4  12  3 

centre,  as  is  T 


are  zero  at  the 


Equation  (23)  reduces  lo:- 
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1  I 


“l2  12 


'66 


^66  ^66  ^^3. 


(28) 


"^2  . 

and  since  and  - —  are  proportional  to  the  applied  acceleration 

d  ^ 

can  write:- 


,  we 


s ,  =  k , ; 


(29) 


ou  . 


=  k^;-  , 


(30) 


The  relative  frequency  variation  is 


then :  - 


^^64'^! 


^  2S6'^2^' 


(31) 


This  follows  a  linear  law  and  shows  that  its  origin  is  due  to  the  effects 
of  shear.  It  also  explains  why  a  change  in  the  direction  of  the  acceleration  F 
produces  a  sign  change  in  the  frequency  variation^. 

Chapter  2  Influence  of  a  constant  acceleration  on  the  resonator 

The  acceleration  applied  to  the  quartz  crystal  produces  internal  forces  and 
hence  surface  forces  through  the  intermediary  of  reactions  caused  at  the  crystal 
supports^’^  by  the  fixing  points.  These  forces  lead  to  static  deformations  of 
the  quartz  which  modify  its  elastic  properties  and  also  its  dimensions.  This 
causes  a  variation  in  the  natural  frequency  of  the  resonator. 

2 . I  Measurement  principle 

When  the  quartz  crystal  is  at  its  resonant  frequency,  the  phase  difference 
between  the  current  in  the  resonator  and  the  voltage  applied  to  its  terminals  is 
zero  (provided  that  the  effect  of  the  parallel  capacitance  and  the  parasitic 

capacitances  are  negligible).  Any  frequency  change  leads  to  a  phase  change,  and 
vice  versa.  The  principle  of  the  method  of  measurement  is  that  of  detecting  this 
phase  change,  using  a  phasemeter  which  gives  out  a  dc  voltage  proportional  to 
the  phase.  This  voltage  is  used  to  servo  the  synthesiser,  externally,  by  means  of 
its  interpolation  oscillator.  The  quartz  crystal  is  excited  at  its  resonant 


12 


frequency  by  the  synthesiser.  However,  when  the  natural  frequency  of  the  quartz 
is  altered  the  corresponding  phase  change,  which  is  detected,  allows  the 

synthesiser  to  be  servoed  to  a  new  value  for  ttie  quartz  resonant  frequency.  The 

frequency  variations  can  be  measured  directly  by  using  the  synthesiser  interpola¬ 
tion  oscillator.  The  phase  changes  can  also  be  recorded  graphically  on  a  graph- 
plotter. 

2 . 1 1  Experimental  equipment 

A  centrifuge  was  used  to  provide  the  continuous  acceleration,  allowing  up 
to  100  g  (g  is  the  unit  of  acceleration)  to  be  applied  to  the  crystal.  Because 
of  the  mechanical  strengths  of  the  quartz  crystal  supports  the  tests  carried  out 
were  limited  to  50  g  in  most  cases. 

The  test  equipment  is  shown  schematically  in  Fig  1.  The  quartz  crystal 
was  excited  by  means  of  the  synthesiser  too  close  to  its  resonant  frequency.  The 

level  could  be  adjusted  with  an  attenuator.  The  power  amplifier  had  a  very  low 

output  impedance,  of  1  H,  and  provided  a  voltage  which  was  practically  constant 
at  the  quartz  crystal  terminals  regardless  of  the  current.  The  resistance 
R  allowed  the  current  in  the  quartz  crystal  to  be  measured;  its  value  was  about 
50  12.  At  the  resonant  frequency,  the  designed  resistance  of  the  quartz  was  of 
the  order  of  several  tens  of  ohms.  The  crystal  was  placed  in  a  chamber  thermo¬ 
statically-controlled  at  the  temperature  of  its  inversion  point.  The  temperature 
stability  of  this  chamber  was  within  0.01°C  per  hour. 

This  experimental  equipment  also  allows  the  parameters  of  the  quartz 
crystal  to  be  measured.  At  resonance,  the  phase  difference  between  A  and  B  is 
zero,  the  voltage  V  is  a  maximum  and  the  impedance  of  the  quartz  crystal  is 

D 

equal  to  the  resistance  of  its  equivalent  circuit  (Fig  2). 


Measurement  of  the  bandwidth  at  3  dB  gives  the  values  of  the  quadrilateral 
frequencies  f^  and  f^  .  The  selectivity  coefficient,  Q  ,  can  be  written 
as :  - 


f 


R|  +  R 


Now 


L  0) 
10 


and  hence 
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K  +  R 


R|  -  K 

27(i7^ 


2 


Cy  is  tlic  intoroleclrodo  capacity  ami  can  bo  moasurod  directly  wlion  tlio 
quartz  is  at  a  distance  from  its  resonant  frequency.  The  parallel  frequency 


can  also  be  used  to  calculate 


('  =  r  _ 2 _ 

0  1  2(f  -  f  ) 

p  s 


where  f  is  the  series  frequency  (f  —  f  ). 
s  ’Sr 

When  tlie  acceleration  is  applied,  tlie  interelectrode  capacity  can 

vary  with  crystal  deformation.  I'his  variation  can  be  related  to  that  of  tlie 
quartz  frequency. 

If  we  let  be  the  capacity  lor  which  the  series  frequency  is:- 


lii  A.  C 


^0  ”  ^0  *  ^  capacity  for  wliich  we  liave:- 


2T,/rTr  '  "  2L, 
(/  1  1  \ 


lew  value  of  the  series  frequency  for  a  given  acceleration. 

“  *-()  ~  i®  related  to  -y  =  by  tlie  equation: 


‘■’•—2 -  ‘'o 
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I 

> 
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To  satisfy  the  experimental  variations  of  f^  (some  tens  of  Hz),  the 
corresponding  variations  of  must  be  of  the  order  of  several  hundreds  or 

thousands  of  picofarads.  Now  is  of  the  order  of  several  picofarads.  Con¬ 
sequently,  we  can  exclude  the  effect  of  as  far  as  the  variation  in  frequency 

of  a  quartz  crystal  subjected  to  a  continuous  acceleration  is  concerned. 

Servo  principle 


This  equipment  (Fig  3)  is  designed  to  achieve  the  quartz  resonant  frequency 
very  rapidly  and  to  follow  the  changes  in  it.  For  that,  we  used  the  synthesiser 
interpolation  oscillator  which  had  a  linear  response  as  a  function  of  the  applied 
voltage. 

The  modulation  slope  of  the  synthesised  frequency  is  a  function  of  the 
modulation  pulsing.  We  chose  a  frequency  range  cf  ±10  Hz  for  a  voltage  of  ±5  V, 
and  hence  for  Af  =  Sv  we  have  S  =  2  Hz/v.  The  time-constant  of  the  inter¬ 
polation  oscillator  is  t  =  l/2iTf  ,  where  f  is  the  cut-off  frequency  relative 

c  c 

to  the  synthesiser  modulation  frequency  defined  at  3  dB. 

f  =  100  kHz,  from  which  1=2  ps. 

c 

The  phasemeter  gives  out  a  dc  voltage  v  proportional  to  the  phase  $ 
measured  at  the  resonator  terminals.  To  measure  the  phasemeter  time-constant  we 
applied  a  reference  signal  (5  MHz  from  an  atomic  clock)  to  one  of  the  channels 
and  to  the  other  channel  a  frequency  signal,  initially  the  same  but  capable  of 
being  varied  slowly.  The  variation  in  the  output  voltage,  proportional  to  the 
phase,  is  shown  in  Fig  4. 

The  operation  of  the  phasemeter  can  be  represented  by  the  following 

g 

differential  equation  :- 

R'C'v'  +  v  =  Kef)  (39) 


where  R'C'  represents  the  time-constant  of  the  apparatus  and  v  is  the  signal, 
a  function  of  d)  . 


The  frequency  for  which  the  phasemeter  output  signal  is  attenuated  by  3  dB 
is  equal  to  about  500  Hz.  The  different  measurement  ranges  for  the  phase  cause 
the  coefficient  K  to  vary.  The  range  used,  +6°  for  v  =  ±  0.5  V  is  the  most 
sensitive  and  determines  the  value  of  K  . 

K  = 


V/rd  . 
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The  servo  sysLem  evidently  does  not  make  an  automatic  search  of  tlie  quartz 
resonant  frequency.  This  must  he  searched  for  by  a  slow,  manual  sweep,  whicli  is 
coupled  to  the  apparatus  when  the  servo-field  is  entered. 

At  this  stage,  the  transitory  slates  are  only  attributable  to  those  of 
quartz. 

2.111  Experimental  results 

The  accelerations  were  applied  systematically  for  all  positions  of  the 
quartz  crystal  (along  the  ox'y'z'  axes)  with  respect  to  the  direction  of  the 
acceleration.  The  tests  were  carried  out  with  different  mountings  having  2,  3, 

4  or  5  supports  for  the  crystal.  Fig  5  represents  the  resonator  in  its 

crystal lograpliic  axes  ox'y'z'.  We  indicate  the  slightly  convex  face  of  the 

crystal,  which  is  a  plano-convex  lens  of  15  nun  diameter  and  of  thicknesses 

I  nan  and  1.65  mm.  All  the  resonators  used  were  cut  in  the  AT  form,  which  is 

the  one  most  used  for  the  production  of  ultrastable  oscillators.  Moreover, 

this  cut  has  the  advantage  of  being  only  sliglitly  sensitive  to  thermal  variations. 

The  plane  face  of  the  quartz  crystal  is  in  the  x'oz'  plane  and  the  convex 
face  along  oy'.  The  supports  were  fixed  inmiaterial 1 y  in  the  z'oy'  and  z'ox 
planes  (Fig  6). 

In  Fig  6  the  figures  in  circles  on  each  of  the  arrows  refer  to  the  synthesis 
curves  (Fig  7)  of  the  frequency  variation  obtained  in  each  case.  The  experimental 
curves,  which  follow,  carry  the  same  reference  numbers. 

2.111.1  Resonator  with  two  supports  (No. 01,  No. 5048  ENSCMB) 

The  curves  of  Figs  8,9  and  10  represent  the  frequency  variations  when  the 
acceleration  applied  was  in  24  different  positions  relative  to  the  crystal. 

With  respect  to  each  of  the  two  senses  defined  by  one  axis,  the  quartz  crystal 
occupies  four  positions,  displaced  by  90*'^,  in  a  plane  perpendicular  to  this  axis. 

It  is  established  that  the  frequency  variations  always  have  the  same  sign 
for  a  given  direction  and  sense  of  the  acceleration. 

Wlien  the  sense  of  the  acceleration  is  changed,  whilst  keeping  its  direction 
the  same  (cases  1  and  2,  or  3  and  4,  or  5  and  6  of  Fig  6)  a  change  in  the  sense 
of  the  frequency  variation  is  noticed  and  in  all  cases  it  can  bo  stated,  defi¬ 
nitely,  that  this  variation  is  directly  proportional  to  the  applied  acceleration. 
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2. 111. 2  Resonator  with  three  supports 

This  quartz  crystal,  in  a  metal  box,  and  of  industrial  manufacture,  shows 
(Fig  II)  curves  of  the  variations  in  frequency  of  which  the  main  characteristic 
is  their  high  degree  of  linearity. 

2. 111. 3  Resonator  with  four  supports 

This  quartz  crystal.  No. 5057,  specially  made  and  held  by  four  supports, 
shows  a  behaviour  somewhat  different  as  far  as  the  relative  values  of  frequency 
variations  are  concerned  (Fig  12).  This  effect  is  due  to  the  increase  in  the 
number  of  quartz  crystal  attachment  points. 

Fig  13  shows  the  results  from  an  industrial  quartz  crystal,  supported  at 
four  points,  which  we  subjected  to  an  acceleration  of  100  g. 

It  would  appear  that  the  behaviour  was  similar  to  that  of  the  previous 
quartz  crystal  up  to  50  g,  but  differed  above  that  value.  The  linearity  dis¬ 
appeared  and  the  shape  of  the  curve  between  50  and  100  g  can  be  attributed  to 
nonlinear  effects.  The  different  points  correspond  to  increasing  excitations  of 
the  resonator.  There  would  appear  to  be  no  behaviour  peculiar  to  the  large 
excitations. 

2.111.4  Resonator  with  five  supports.  No. 5065  (Fig  14) 

This  resonator  was  constructed  specially  at  the  ENSCMB  piezoelectric 
laboratory  in  order  to  study  the  action  of  a  large  number  of  supports  as  well 
as  the  influence  of  the  positions  at  which  these  supports  are  attached.  It 
would  appear  that  the  frequency  inversion  disappears  for  the  oy'  axis  and  that 
the  relative  values  are  changed  with  respect  to  the  first  results.  All  the 
available  evidence  shows  that  the  anchorage  points  of  the  supports  play  a  large 
part  in  the  quartz  crystal  behaviour.  It  seems,  from  this  example,  that  a  care¬ 
ful  choice  in  their  positioning  can  give  compensation  for  certain  effects. 

Remarks 

(1)  In  the  case  of  resonators  with  two  supports  -  the  type  least  disturbed  by 
the  presence  of  fixing  points  -  it  was  noticeable  that  acceleration  along  the 
ox  axis  produced  the  smallest  frequency  variation.  Now  ox  (or  ox')  is  the 
only  axis  which  keeps  to  the  original  crystal  symmetry.  It  seems,  then,  that 
the  disturbance  effect  is  the  least  sensitive  when  it  is  applied  along  a  parti¬ 
cular  axis  of  the  quartz  crystal. 

(2)  In  the  case  of  a  quartz  crystal  with  two  supports,  we  take  a  mean  value  of 
the  frequency  variations.  The  variation  is  reasonably  linear,  and  we  can  write:- 
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Af  =  kl 


and  Lhe  values  of  k  for  curves  1  to  6  of  Figs  8,  9  and  10  are:- 


kj  =  -  3  H  10  "  Hz/g 
k^  =  A. 4  X  I0~^  Hz/g 
k^  =  -  1.2  X  10~^  Hz/g 

k,  =  1,2  X  10~^  Hz/g 

k^  =  -  1.6  X  10~^  Hz/g 

k^  =  2  X  lo"^  Hz/g  . 

The  maximum  value  of  the  coefficient  k  was  obtained  for  a  quartz  crystal 
with  five  supports  and  was  0.04  Hz/g. 

2. III. 5  Variation  of  the  anti-resonance  frequency 

In  order  to  provide  a  better  definition  of  the  quartz  crystal  behaviour, 
the  same  measurements  were  carried  out  on  the  parallel  frequency.  However, 
this  is  not  very  pronounced  for  a  quartz  crystal  with  a  small  interelectrode 
capacity  .  By  placing  a  capacity  of  about  10  picofarads  in  parallel  with 

the  quartz  crystal,  the  an t  i-rcsonance  frequency  is  brought  close  to  the  resonant 
frequency  and,  in  particular,  the  minimum  current  in  the  quartz  crystal  is  then 
much  more  easily  displaceable. 

The  variations  in  the  anti-resonance  frequency  maintain  the  same  character¬ 
istics  as  those  found  previously  for  resonance. 

Figs  8,  9,  10  and  1 1  show  these  variations  for  the  quartz  crystals 
considered . 

By  using  the  relationships  between  the  resonance  (wq)  and  anti-resonance 
(u)  )  frequencies  it  is  possible  to  check  that  these  frequencies  are  identical. 

3 


hVi 


we  have 
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(D  «  U)  (1  +  - — 

a  0  ^  2CJ  • 


Let  be  the  new  value  of  during  acceleration,  is  of  the  form 

u)*(l  +  6^),  where  u)*  is  the  constant  value  of  the  anti-resonance  frequency. 


Simi lar ly : - 


l!  =  c*ii  +  n) 


c'  =  c*u  .  ,) 


and  equation  (40)  is  equivalent  to:- 


.*(1  ^  V 


C* 

+  —L  '  h 

2C*  1  .  , 


and  hence:- 


“a  "  "O  I  '  "  ^  "  n  -  >)  -  ^ 


U) '  ^  U)*  (  I  + 

a  0  V 


u.*S*  # 

a  a 


from  which 


2C*J  ^  “"O  2C*  ^  ^  "’O^O  ^  ^  h  -  ^) 


c* 

til*  6  +  — L  r 

o|^0  2C*|_ 


(n  -  >)  +  6^(1  +  n  -  >) 


m*  [  0  *  2C*  [' 


(n  -  >)  (I  +  iS^j)  -»■  6^^ 


Now  6q  is  small  with  respect  to  1  and  n  and  >  are  very  small  with  respect 

and  also  C*/2C*  is  of  the  order  of  10  .  Consequently  it  is  justifiable 
to  write:- 


w*  (1 


It  can  be  seen  tliai  A  is  practically  equal  to  the  variation  6^^  of  tliv 
lesonant  frequency  and  so  this  calculation  confirms  tlie  experimental  results. 

All  the  variations  in  the  resonant  and  ant  i-resonaiue  frequencies  show  that 
the  effect  of  the  acceleration  on  the  quartz  crystal  is  to  cause  a  translation  of 
the  'current  in  the  (juartz  crystal'  curve  as  well  as  of  the  phase  curve 
(fig  15). 

d.lll.b  Effect  of  hysteresis  ol  the  lre([uency  variation 

When  the  acceleration  decreases  the  frequency  does  not  go  back  to  the  sanii 

values  which  it  had  when  the  increasing  acceleration  was  applied.  No  complete 

explanation  for  this  is  given  here,  but  it  is  comparable  with  the  phenomena  slKn^ 

. .  .  .  S)  . 

when  a  field  is  applied  .  It  is  probably  due  to  a  progressive  disappearance  ot 
the  stresses  caused  by  the  memory  effect  of  the  crystalline  structure.  When  tlu 
acceleration  is  zero  once  more,  the  difference  between  the  initial  and  final 
trequencies  can  reach  0.1  Hz  and  it  is  jiossible  for  the  variation  to  change  sign 
during  the  course  of  the  deceleration  (I'ig  lb).  It  often  takes  several  tens  of 
seconds  for  the  resonator  to  return  to  within  one  hundredth  of  1  Hz  of  its  initi 
t  requency . 

d.111.7  Variation  in  the  designed  resistance 

In  order  to  determine  the  equivalent  resistance  of  the  resonator  it  is 
sufficient  to  measure  the  voltage  at  the  terminals,  when  the  resonator  is 
excited  at  its  resonant  frequency.  We  then  have: 


V.  -  V 

P  .  A  B  p 

‘'i  “  — V  '' 

B 


The  voltages  and  are  measured  by  means  of  a  phasemeter.  If 

and  AVj^  are  the  variations  of  these  voltages,  respectively,  then  AR^  , 
the  variation  in  the  designeil  resistance,  is  such  that:- 


AR,  0  if 


and  vice  versa. 


Moasuromoius  show  thal  i«  all  eases  ^  s  ^  . 

j  i  ,  .  Va  '  from  this  it  can  be 

euooa.ha,  t.„.  K,  a."  tl,.  a"cco i „„  i 

:t::: . \ ^  -« va,„.. ... . .. 

U.5  1.  (tor  a  resistance  of  I00  i!)  at  50  g. 

,,  „„  . . 

to  U,.  ,. 

••  -.•  ...otto,,  ,„o  o,.e-.i„.,.to,.ot. 

.«  o.  tH.o  „o,to„  t.  to  o,,toi,.  t„o  ottoot  o,  t„o  ootvot,,..  oot.eo. 
large  accelerations. 

We  use  the  fuiid.iraent al  equations  (Refs  10  II  in 

VMis.iu,  II,  12)  anwngst  which  are:- 

the  equation  for  dynamic  equilibrium 

") 

^  u . 

1 

P  ~  =  d  1 V  T .  +  K 

3t"  J  i 

.<-.'1— nts.  otroooeo  o„.  t,.. 

the  force-deformation  equation: 


AV,  AV., 


T.  »  CV.S.  +  -  ( '  s  q  c  _  ..  '"’i 

‘  ij  j  h  ijke^j^k  C  ‘’nu\,  ^  '’ij  (43 

"litre  and  represent  the  deformations  and  electric  fields. 

.,  /ij’  ^ijk<’  ‘mi  '‘ij  Uie  Jnd  and  4th  order  elastic  coefficients 

‘  Piezoelectric  coefficients  and  the  damping  coefficients,  respectively. 

the  electrical  displacement  1) 


e  .S.  +  t-  E 
iiJ  J  mn  m 


^mn  eoefticient  of  electrical  polarisation. 

We  used  the  mathematical  models  developed  bv  1  I  r . ,  • 

bv  him  to  the  nt  •  •  ^  and  applied 

the  determination  of  the  influence  of  a  continuous  field  on  the 

resonant  trequency  of  a  quartz  crystal''^. 


In  t  lu“  caso  ol  a  plate  iiil  along  AT  ,  with  thickness  e  ,  ol  infinite 

all  al  li  inii'iisions ,  vibrating  in  shear  across  the  thickness,  the  deformation 

.  iiii  es  to  ,  the  stress  to  1'^^  and  i  hi'  lield  is  aiijilied  along  oy  .  .Ve 

hall  neglect  the  damping  r.  .  whiiTi  does  not  cause  anv  non  1 i near i t i es . 

ij 

hqiiations  t4J),  1)  and  (.4a)  can  he  written:- 


dl 


dy 


T  =  C'^S 
(>  b  b 


^  —  C  S  -  e K 
b  bbt>b  b  Jb  2 


2  b  b 
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N  .in  obiain  tho  propagation  oquation:- 


0  u 
> 

3t" 


bb  4  2  bbbb 

3v 


+ 


pl' 


hiih  are  associated  the  two  tollowing  limiting  conditions:- 


u  =  0  for  y  =  0 

T  =  0  for  y  =  * 

o 


i 


'Ol 


If  we  Lake  a  solution  of  the  form:- 

u  =  *■  A(,y)  cos  u't  i.-iO'i 

1),  =  11^^  +  My  cos  lot  .>'■ 

l.f 

184/  an-  pnt  this  into  ecpiation  (44)  we  obtain  the  following  equations:- 


I;  ; 
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-  1  ^6666  ^  ^  ^0  f ^  ^  ^  ^ 

2  C°.  Idv^  W  /  2  \dy)  *  2  dy  dy 


^2.  2 

2  2  * 

dy  c 


1  6666 

2  C° 

66 


h.,2  dy  dy  .,2  I  dy 


1 

^  ^  dy^ 


:5::  .  (55) 


In  each  of  these  two  equations  the  second  term  is  small  and  allows  solutions 
to  be  obtained  by  successive  approximations.  Thus  we  obtain:- 

Un  =  (±ey  -  y^)  (56) 


which  is  the  general  solution  for  equation  (54)  without  the  second  term,  and  is 
accurate  enough  for  this  calculation. 

In  the  same  way  we  have:- 


A(y)  =  a  sin  —  y 
c 


and  if  we  put  equations  (56)  and  (57)  into  equation  (55)  we  obtain:- 


2  2 

+  = 

,2  2  * 

dy  c 
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\^66/ 
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~  r~D“  "  2y)  sin  -  y 

VS6/ 


3  3  to  .  (0 


r-  a  — r-  sin  —  y  cos 
4  4  c 


whose  general  solution  is:- 


A  =  a  sin  ^  y  -  ^ 

c  ^  2  D 

^66 


666  /  .  toy  loy  2  (DV 

—  (  “,y  sin  +  a^y  cos  +  a^y^  sin  ^  +  a 

66  ^ 

^  3  ojy  .  .  UJY  \ 

+  a  y  cos  — ^  +  a  sin  e  — ^ ) 
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where 


(60-1) 


(60-2) 


(60-3) 


(60-4) 


(60-5) 
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The  application  of  the  second  limiting  condition  leads  to  the  following 
equation:- 


'26 


66  6  6  ^6666^6  ^2 


=  0 


for  y  =  ±  ie  (61) 


and  from  this  we  have:- 


C  _  ^  ^  ^1 


If  we  put 


'26 
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^6  =  ^^‘^2 


(63) 


and 


'  4 

f 


we  can  obtain,  by  identification:- 


-  I 


3  3  2 


aD  -  —  eD  =  — ~  I*  —\  = 

0  Wo  2  0”0  dy  2 ) 


where 


^  6  ^6666“ 


By  using  equations  (59)  and  (60-1)  to  (60-6)  equation  (64)  can  be  expressed 
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for  y  =  ±  -e  (67) 
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The  conditions  which  affect  the  field  must  then  be  applied; 


D  e 
2  26  3u 


2  £22  ^22 


By  integrating  over  the  crystal  thickness  we  obtain:- 


After  identification,  two  new  equations  appear: 
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and 
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By  taking  into  account  Che  orders  of  magnitude  of  the  different  coefficients, 
equations  (67)  and  (72)  reduce  to;- 
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(74) 


Since  we  are  not  interested  in  the  absence  of  isochronism,  that  is  to  say 

3 

in  the  influence  of  the  level  of  oscillation  on  the  frequency,  the  terms  in  a 
can  be  neglected.  The  dependence  of  the  frequency  on  the  applied  acceleration 
will  be  obtained  by  letting  tend  towards  zero.  We  then  obtain:- 
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where  6^^  represents  the  relative  error  between  the  excitation  frequency  and  the 
natural  frequency  of  the  resonator,  and  n  the  range  used  in  the  vibration 
(here  n  =  5) . 

is  clearly  a  quadratic  function  of  the  acceleration. 

Chapter  3  Influence  of  a  sinusoidal  vibration  on  the  resonator 

The  problem  of  vibrations  is  different  from  the  previous  one  in  the  sense 
that  it  introduces  dynamic  deformations  of  the  crystal  and  hence  a  low-frequency 
modulation  of  the  resonator  natural  frequency. 

3 . I  Measurement  principle 

The  method  used  consists  of  detecting  the  frequency  variations  by  means  of 
the  corresponding  phase  fluctuations  in  the  current  through  the  resonator.  For 
that,  a  spectral  analysis  is  carried  out  on  the  output  signal  from  the  phasemeter. 
Since  the  phase  fluctuations  are  random,  it  is  possible  to  extract  from  this  noise 
the  harmonic  components  in  direct  relation  to  the  sinusoidal  disturbance. 

3.11  Experimental  equipment 

The  experimental  equipment  is  almost  the  same  as  that  used  for  steady 
accelerations,  but  the  manipulations  were  carried  out  without  the  frequency  servo- 
loop.  In  this  case  we  used  the  voltage  output  from  the  phasemeter  connected  across 
the  terminals  of  the  quartz  crystal.  This  voltage  was  recorded  on  an  FM  tape 
recorder  then  analysed  on  a  programmable  analyser  controlled  by  a  computer.  The 
accuracy  of  the  analysis  depended  upon  the  freqaency  range  chosen.  To  give  an  idea 
of  the  order  of  magnitude,  the  incremental  accuracy  was  1  Hz  for  a  500  Hz  range  of 
analysis.  The  analysis  dynamic  was  50  db  and  each  of  the  spectrum  graphs  was 
accompanied  by  a  table  showing  the  frequency  and  amplitude  of  each  line  with 
respect  to  a  reference  signal. 

3.111  Method  of  analysis 

It  is  possible  to  analyse  a  certain  frequency  range  in  real  time  thanks  to 
compression  of  the  signal  in  time  and  which  also  allows  the  time  taken  for 
analysis  to  be  'accelerated'. 

The  analysis  is  carried  out  in  an  analogue  matter  but  on  the  other  hand  the 
signal  acceleration  is  achieved  by  the  use  of  numerical  techniques  (shift 
register,  acoustic  delay  line). 

To  explain  the  principle  of  this  analyser,  the  real  acceleration  mode  is 
replaced  by  a  sjmibolic  acceleration  mode.  Suppose  that  a  portion  of  a  signal 
x(t)  of  duration  At  is  recorded  with  a  writing  speed  v  on  a  magnetic 


tape- loop.  Tlu'  signal  can  tlicn  bo  road-ol  I  by  moans  ol  a  roail-o\il  lioad  monnt»'il 
on  an  arm  and  turning  at  a  spood  V  »  nv  (a  1).  Tlu'  signal  x'(t)  road 
during  t  bo  first  rovolution  ol  t  bo  roaii-out  boad  will  bo  t  bo  accoloratoii  signal 
and :  - 


x‘  (t) 
a 

"  x(ut) 

who  re 

u  is  the  acceleration  ratio. 

If  we  put 

X(v) 

-  TFx(t) 

where 

TK  is  the  Fourier  transform 

t  hen 

we  h.ive:- 

x‘  (v) 
a 

“  lTx\t) 
a 

A  t  /  u 

r 

=  j  xlut)e 

6 

«=  —  Xlv/u) 

Tbus  Lbo  spoctral  sproad  ol  a  poriod  ol  an  acooloratod  signal  and  that  ol 

tbo  initial  signal  arc  in  tbo  ratio  a  .  In  ordor  to  I'btain  tlio  spootrnm  I'l  i  I  <■ 

signal  x(t)  witbin  tbo  f ro<iuoncy  rango  111, Nil)  by  moans  ol  N  oqnidistant 

tiltoring  points,  and  with  a  resolution  11  =  I/At  ,  it  is  sniticioni  to  analys. 

X  (v)  with  a  filter  of  resolution  ull  . 
a 

Again,  tbo  filter  and  tbo  associated  detection  system  have  a  response  t i  im 

^,1  —L  <E  AL  wbicb,  for  signal  analysis,  allows  tbo  use  ol  a  sliding  liltoi  !■ 
ull  u 

whose  central  fro<iuoncy  v^.  takes,  sviccessively,  the  values  ull,  /ull  ...  Null  at 

,  ,  At  (N  -  1)  . 

the  lime  instants  ol  1),  —  ...  -  At  . 

u  u 

/\  t. 

The  analysis  lime  T  -  N  —  is  equal  to  the  sign.al  duration  in  N  «  u  . 

ll 

and  the  system  operates  as  if  the  signal  xlt)  ,  of  dural itin  At  ,  bad  been 
ana ly soil  by  N  filter  points.  In  reality,  the  signal  xU)  is  api'lied 
continuously  to  the  analyser  in|<ut,  which  gives  out  an  elementary  spectrum 
|x*'lv.)|  ,  ij  ”  I  to  N  ),  the  spectrum  obtained  during  the  kth  antilysis  seipi.  <oe. 
This  spectrum  passes  through  a  qutidratic  circuit  and  then  its  mean  value  ovei 
I  reipiency  range  is  found. 


i.r 
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:(l) 


c\v.) 


S(v  .  ) 
J 


(S  (v.)  spectral  luuctiuu  is  thus  related  to  |x  (v . )  1 
by 


0  (v 


(77) 


k=l 


S(v)  ,  tlie  power  spectral  density  associated  with  x(t)  is  such  that 

S(v)  -  XE  l]  where  X  is  a  scaling  constant  peculiar  to  the  apparatus. 

Hence  G*^(Vj)  is  a  correct  estimate  of  the  power  spectral  density  S(v)  . 

i. IV  Experimental  results 

We  have  applied  to  tlic  resonators  sinusoidal  vibrations  within  a  frequency 
range  of  10  to  5000  Hz  and  from  several  g  to  40  g.  The  quartz  crystals  were 
placed  in  a  specially-constructed,  thenmjstatically-controlled  chamber  to  undergo 
the  vibrations  (Eig  17). 

The  resonators  used  could  form  different  sources,  thus  enabling  the  testing 
of  quartz  crystals  having  two,  three,  four  or  five  support  points.  For  low- 
frequency  external  vibrations  the  effect  of  the  supports  is  due  to  their  moments 
of  inertia.  Only  those  resonators  manufactured  by  the  ENSCMB  laboratory  were 
subjected  to  the  highest  accelerations.  The  quartz  crystals  with  two  supports 
formed  part  of  the  normal  laboratory  manufacture  but  those  with  four  and  five 
supports  were  manufactured  especially  for  this  study. 

5 . 1 V .  1  Manipulation  of  the  frequency  servo  with  closed  loop 

Although  frequency  servoing  was  not  necessary,  a  series  of  tests  was 
carried  out  under  these  conditions.  The  vibration  frequencies  were  chosen  spon¬ 
taneously  at  any  suitable  values,  but  different  from  those  of  the  sector  and  of 
its  harmonics.  Fig  18  shows  the  fluctuations  in  the  natural  frequency  of  quartz, 

S>  ,  when  the  quartz  crystal  does  not  undergo  any  disturbance  and  Figs  19  to  21 
0 

show  the  fluctuations  for  exciting  frequencies  and  accelerations  of:- 


12  Hz,  3  g  Fig  19 

bO  Hz,  3  g  Fig  20 

87  Hz ,  3  g  Fig  2  1 . 
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The  appearance  of  harmonics  of  the  vibration  frequency  can  be  seen,  which 
can  only  be  explained  as  due  to  the  servo  loop.  With  the  loop  open,  the  harmonic 
peaks  disappear.  In  fact,  when  the  loop  is  closed,  the  synthesiser  is  guided  by 
a  voltage  whose  harmonic  content  (fluctuations  in  the  resonator  frequency)  is 
regenerated  within  the  synthesiser.  Each  spectrum  shows  the  presence  of  a  large 
peak  of  fixed  frequency  (particularly  Fig  18).  This  arises  from  the  mixing 
between  the  natural  frequency  of  the  quartz  crystal,  shown  by  the  synthesiser  and 
the  5  MHz  calibration  frequency  which  guides  the  synthesiser. 

This  remark  allows  us  to  deduce  the  relative  amplitude  of  the  disturbance 
peak  with  respect  to  the  carrier. 


If  we  consider  the  scheme  above, 
where  Wq  is  the  quartz  crystal  resonant  angular  frequency 
0)^  is  the  reference  angular  frequency  =  2it  x  5  MHz. 


At 

A 

we  have 

=  ^1 

sin  Wpt  +  Vj 

sin  0)  t 
r 

and 

at 

B 

have 

"b 

=  ^^2 

sin  Wq!  + 

sin  0)  t  • 
r 

The  resonators  have  a  quality  coefficient  of  1.5  x  10  to  2  x  10  and  the 
5MHz  component  is  practically  nonexistent  at  the  quartz  crystal  output,  since 
this  plays  the  role  of  a  filter,  and  so  can  be  neglected. 

The  phasemeter  gives  an  output  voltage  comparable  with  that  from  a  4-diode 
ring  demodulator. 


V  =  V  V  =  ■  '*  ^  cos  Amt  +  g(2u  wr  +  io„) 

s  A  B  2  0  0 

where  Aw  =  “j.  “  the  difference  between  the  angular  frequencies; 

g(2w^wr  +  Wp^)  represents  the  terms  with  a  frequency  of  10  MHz. 

0  0  Vj  -  V2 

Consequently,  -rrpry -  represents  the  ratio  between  the  voltage  at  the 

12 

2 

quartz  crystal  terminals  and  the  parasitic  peak. 


Tho  values  currently  in  use  are:- 


V|  =  JO  niV 
V,  =  10  mV 

Vj  =  500  iiV 

wliicli  ^ive  a  ratii'  ol  12  liB. 

Since  Llie  analysis  ^ives  Llie  error,  in  ilB,  between  t  lie  parasitic  Ireciuency 
peak  and  that  ot  tlie  disturbing;  Ireciuency,  it  is  simple  to  obtain  tbe  ratio  between 
tbe  amplitude  ot  the  ciuartz  crystal  resonant  lrec|uency  and  that  of  tbe  external, 
disturbing;  vibration. 

As  in  the  case-  of  the  steady  accelerations,  three  particular  cases  will  be 
cons i de  red . 

) .  1 V .  J  Direction  of  vibration  parallel  to  ox 

this  arrangement  is  such  the  natural  vibration  of  the  crystal  and  the 
disturbing  vibrations  are  in  the  same  direction.  figs  Id  to  J1  show  that  the 
ettect  ot  the  vibration  diminishes  as  the  ireciuency  increases  and  that  the 
icceleration  remains  constant.  In  the  case  where  the  ireciuency  is  fixed  and  the 
acceleration  varies,  the  amplitude  of  the  peak  increases  with  the  acceleration, 
as  Figs  JJ  to  J4  show,  where  the  acceleration  has  the  values  of  J  g,  4  g  and  8  g, 
respectively,  at  a  Ireciuency  of  30  Hz. 

I'liese  results  are  the  same  whether  the  external  vibration  is  directed 
long  the  +ox'  axis  or  along  the  -ox'  axis  of  the  crystal. 

1 . 1 V . J  Direction  of  vibration  parallel  to  oy' 

In  this  configuration  the  two  faces  of  the  crystal  are  perpendicular  to  the 
vibration.  Again,  the  position  of  one  particular  lace  with  respect  to  the 
disturbance  is  of  no  importance.  It  appears,  however,  that  in  this  position  t he 
resonator  is  much  less  disturbed  than  in  the  previous  case.  Moreover,  fc'r 
increasing  acceleration  at  a  fixed  freciuency,  the  disturbetnee  is  practically 
unchanged  (Figs  25  and  2b).  However,  it  decreases  when  the  freciuency  increases. 

3 .  I V . 4  Direction  of  vibration  parallel  to  oz* 

This  position  gives  results  which  are  in  all  ways  comparable  to  those  for 
the  first  case  and  the  conclusions  are  the  same. 


51 


Besides  the  phase  inlonnation,  the  phasemeter  gives  a  measurement  of  the 
signal  amplitude.  It  provides  a  dc  voltage,  proportional  to  the  ^unplitude  of 
the  observed  signal,  and  whose  t luctuations  are  represented  by  those  of  the  dc 
voltage.  We  have  proceeded  with  the  analysis  of  this  voltage  without  being  able 
to  detect  the  presence  of  the  disturbing  peak,  which  lias  led  us  to  think  that  the 
phenomenon  seen  for  the  whole  of  the  quartz  crystals  is  a  frequency  iiKidulalion 
phenomenon  without  amplitude  modulation.  Referring  to  Kig  15,  we  see  that  the 
external  vibration  phenomenon  operates  on  the  phase  curve  and  produces  a  dis¬ 
placement  of  the  external  disturbance  trequency.  This  confirms  the  frequency 
modulation  phenomenon. 

Kig  27  summarises  all  the  measurements  carried  out.  Kach  curve  represents 
the  ratio  Aii/A^  in  dB,  that  is  to  say  the  ratio  of  the  amplitude  ol  the  dis¬ 
turbing  peak,  of  frequency  si/2ii  ,  to  the  amplitude  of  vibration  ol  the  quartz 
crystal  at  its  resonant  freciuency.  The  el  feet  of  the  disturbance  increases  as 
the  frequency  drops  and  as  the  acceleration  increases.  The  positions  for  which 
the  ox'  and  oz'  axes  of  the  quartz  crystal  are  vertical  (designated  by 
vertical  positions)  give  larger  disturbances  than  are  obtained  wlien  the  oy'  axis 
is  vertical  (shown  by  the  horizontal  position  of  the  quartz  crystal). 

The  vertical  position  oz'  appears  to  be  more  sensitive  tlian  the  ox' 
position  for  low  frequencies.  This  is  due  to  the  intluence  of  the  supports, 
which  disappears  above  100  Hz. 

Technological  limitations  in  the  exciter  did  not  allow  more  than  ^  g,  at 
10  Hz,  to  be  applied.  Also,  the  dynamics  of  tlie  analyses  (50  dB)  has  limited 
the  researches  to  disturbance  frequencies  greater  than  500  Hz  and  for  which  the 
effect  is  extremely  small.  However,  we  shall  show  that  they  show  up  clearly  in 
the  oscillator  phase  spectrum. 

We  have  tried  to  obtain  a  mathematical  description  ol  the  curves  relative 
to  the  vertical  position  ox'  by  using  the  method  of  least  squares.  In  an 
orthogonal  reference  they  can  be  written  in  the  form:- 


a.r^  +  b.r 
1  1 


(7B) 


whe  re  a .  ■»-  b  . 

i  1 
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W.  IS  fixed  for  each  curve  and  Fig  28  shows  the  variations  in  a.  and  b. 
astunctions  of  the  vibration  frequency  Q  .  The  slope  of  this  straight  line  is’ 

1  ‘  •  Consequently,  the  variation  of  the  resultant  disturbance  amplitude 

undergone  by  the  quartz  crystal  obeys  a  1/f^  law.  The  same  results  are  true 
for  the  vertical  position  oz'  . 

'  .V  Conclusions 

Whilst  being  subjected  to  a  sinusoidal  vibration,  the  quartz  resonator 
gives  out  a  frequency  which  is  modulated  by  the  vibration  frequency.  The  modula¬ 
tion  increases  as  the  vibration  frequency  drops  and  as  the  acceleration  increases. 


el  tec t. 


Finally,  the  position  of  the  crystal  determines  the  size  of  the  disturbance 


Ch_ap_ter  4 _ Vibration  of  the  (quartz  resonator  supports 

The  behaviour  of  the  quartz  resonator  is  modified  by  its  supporting  elements 
which,  by  their  natural  resonance,  can  induce  considerable  disturbances  in  the 
quartz  crystal. 

We  shall  examine  the  quartz  crystal  with  two  supports,  which  is  the 
traditional  method  of  manufacture,  by  considering,  successively,  the  transverse 
bending  vibrations  and  the  symmetrical  vibrations  of  an  arched  structure. 

•  1  Transverse  bending  vibrations 

The  model  used  is  shown  in  Fig  29.  The  mass  M  represents  the  quartz 
crystal  and  the  support  housing  is  animated  by  a  sinusoidal  vibration'^'"’. 

The  general  equation  for  free  transverse  vibrations  of  a  bar  can  be 
wr i tten;- 


El 


3x  L  3x  J 


-  pS 


at 


i7b) 


where  y(x,t)  defines  the  transverse  displacement  of  the  centre  of  the  section 
K  =  Young's  modulus 

=  the  bending  moment  of  inertia 
=  the  density 

=  the  cross-section  of  the  bar. 


In  the  case  where  E  and  I  are  constants,  the  equation  becomes :- 


11 
18  P 
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ps  -  4 
3t  3x 


The  use  of  normalised  functions  Y(x)  gives  a  solution  of  the  fo 


rm:  - 


y(x,t)  =  Y(x)e 


iS2t 


where 

with 


Y(x) 


a  cosh  ax  +  b  sinh  ax  +  c  cos  ax  +  d  sin  ax 


A 


a 


pS 

El 


W7) 


(78) 

(79) 

(80) 


a,  b,  c  and  d  are  constants  of  integration  defined  by  the  limiting  conditions 
imposed  by  the  methods  of  fixing  the  supports. 


In  our  case,  the  support  housing  was  excited  at  the  pulse-rate  12  by  a 
forced-vibration 


y  =  A  sin  fit 


The  limiting  conditions  are:- 
-  at  the  housing 

v(o,t)  =  A  sin  fit  ->  Y(0)  =  A  (81) 


at  the  free  end 


0  -*■  Y'(0)  =  0  for  fa  +  c  =  A 

<a(b  +  d)  =  0 
|a  #  0 

.  (82) 


Y"(L)  =  0 


(83) 


because  the  flexing  moment  is  zero. 

However,  the  force  T  =  -  EIY'”(x)  is  equal  to  the  inertial  force  due 
to  the  mass 


>T 
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and  hence 


K1Y"'(L)  =  -  MirVd.) 

Y"’(h)  =  -  ^  Y(L)  . 

c.  1 

Now:- 

2  4  4 

NUJ  =  _  Mu  L  ^  M  4 

El  pS  pSL  ra  ^ 

where  m  is  the  mass  of  the  lamina 
and 


Y"'(L) 


M  4 

-  a  LY(L) 

m 


184) 


By  putting  p  =  N/m  and  taking  into  account  equation  (82),  Y(x)  becomes:- 


Y  (x)  -  a^cosh  ax  -  cos  ax)  ■+  b(sinh  ax  -  sin  ax)  +  A  cos  ax  (85) 


and  equation  (84)  gives 

a(sinh  \  —  sin  •)  btcosh  >  -*■  cos  ' )  A  sin  \  =  ..\  lad'i'-‘’h  '  ~  cos  \) 

-t  b(sinh  \-sin  \)  A  cos  \  I 

.  (8b) 

where  A  =  aL  , 

Equation  (83)  becomes:- 

a(cosh  \  -t-  cos  \)  b(sinh  \  +  sin  \)  -  A  cos  A  =  0  (87) 

a  and  b  can  be  calculated  from  equations  (8b)  and  (87):- 

a  ~  ^  I  sinh  A  -t  sin  A  sinh  \  cos  A  cosh  A  1 1  (88) 

b  =  ^  ^hA  cos  A  cosli  A  +  cos  A  sinh  A  +  cosh  A  sin  A1  (8d) 
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where 

anil 


A  =  u\(siu  \  cosh  \  -  cos  \  sinh  \)  -  cos  \  cosh  A  - 


(^0) 


y(x,t)  = 


A  sin  t 


^cos  ax  -  cosh  ax)  (,  1  +  sin  \  sinh  A  +  cos  \  cosh  X 


+  V  cos  X  sinh  \)  +  Xsinh  ax  -  sin  ax) (sinh  A  cos  X 


+  cosh  X  sin  X  +  .’iiX  cos  X  cosh  A) 


+  cos  ax  A  Xsin  X  cosh  \  -  cos  X  sinh  X)  -  cos  \  cosh)  - 


X'^l) 


Wlien  the  imposed  piilse-rotc  h  is  such  that  A  =  0,  the  lamina  begins  to 
resonate.  This  happens  when  h  coincides  with  one  I'f  its  normal  free  mode  pulse- 
rates  . 

We  require  the  first  root  of  the  transcendent  equation  X^O)  ,  in  which 
U  =  34,  because  the  mass  of  the  quartz  crystal,  M  =  0.783  g,  is  supported  by 
two  strips  of  mass  m  =  ll.b  mg. 


We  f i nd 


=  0.5-44 


and  from  X  =  aL  and  a"*  = 
system  begins  to  resonate:- 


we  can  deduce 


the  pulse-rate 


for  which  the 


(82) 


We  have  designed  a  mounting  analogous  to  the  quartz  crystal  support  (Fig  30), 
in  which  the  resonator  is  replaced  by  a  circle  of  duralumin.  Since  their  densities 
are  similar,  the  crystal  dimensions  can  he  maintained.  To  this  circle  is  cemented 
a  miniature  accelerometer  of  very  low  mass  (140  mg)  and  the  whole  is  equal  to  the 
mass  of  a  partial  quartz  crystal  5,  that  is  to  say  0.785  g.  The  accelerometer  is 
a  pietzite  which  gives  a  load  variation  at  the  amplifier  input.  An  automatic  sweep 
of  the  excitation  frequency  is  made  and  a  recording  is  made  of  the  output  voltage 
corresponding  to  the  response  of  the  element  under  study  as  a  function  of  frequency. 
The  system  is  shown,  schematically,  in  Fig  31. 


I 

y 

I 


Jb 


As  a  numerical  example,  let;- 


,\y  =  0.544 

L  =  1 4  X  1 0  ^  m 


K 

1 


1  1  2 
2  X  10  N/m 

10"*^  4 

— r:; —  m 


P 

S 


8,9  X  10^  kg/m 

-7  2 

10  m 


1^  and  are  the  moments  of  inertia  corresponding  to  positions  (a)  and  (b) , 

respectively,  of  Fig  29. 


For  position  (a)  we  have,  by  applying  equation  (92):- 


Oa 

2ti 


529  Hz  . 


Fig  32  shows  the  response  curve  for  such  a  case  and  indicates  the  presence 
of  the  main  resonance  peak  at  281  Hz.  This  is  confirmed  by  experiment.  The 
difference  between  the  two  values  is  easily  explained  by  the  choice  of  an 
analytical  model  which  is  not  an  exact  counterpart  of  the  real  case. 

.  ^‘'b 

For  position  (b),  —  =  33  Hz.  This  is  the  position  which  offers  the 
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smallest  moment  of  inertia  to  flexing.  Referring  to  Fig  27,  we  see  that  for 

30  Hz,  the  disturbance  A, /A,,  has  the  greatest  value.  This  shows  the  direct 

il  0 

action  of  the  support  in  such  a  case. 


The  solution  y(x,t)  also  enables  us  to  find  the  pulse-rate  for  which  the 
displacement  of  the  mass  is  zero,  by  solving  the  equation  y(L,t)  =  0  .  We  shall 
not  do  this,  here,  as  it  is  only  of  limited  interest. 

4.11  Symmetrical  vibrations  of  an  arched  beam 


The  assembly  of  two  supports  and  the  crystal  forms  an  arched  beam.  We  have 
studied  its  symmetrical  vibrations  whilst  neglecting  the  longitudinal  and  trans¬ 
verse  vibrations  of  the  strips.  Also,  the  model  used  is  that  of  an  arched  beam 
with  its  load  distributed  over  its  horizontal  part  (Fig  33). 
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The  method  is  to  1  iiui  the  static  del  orma  t  ions  of  tlie  beams  and  to  introduce 
these  into  tlie  system  energy. 

Consider  0  ,J  (Kig  34a). 

To  find  the  static  detlexion  6  in  J  we  superimpose  a  fictional  force  F 
on  the  distributed  load  4  . 

At  H,  ,  the  flexing  moment  is:- 

Mjl  (x)  =  C  -  fl.  -  x)-’  -  K(L  -  x)  . 

The  deformation  energy  W,  is  given  by:- 


“•  •  Tirrr;  I 


2i: 


C-i,  -  CFI,"  ^  ^  FL^  +  ^ 

3  -4  20 


193) 


Considering  |  ^^8  )4b,  let  C^  be  the  couple  resulting  from  the 

action  of  0  ,J  at  0,  . 


and 


C  =  C  -  4  -  FL 


^1 


1 


-  x) 


and 


“■  =  -v7  i 

0 


'  2F 


.  (4  .  K)[L-(:f=  .  k)  -  ;a| .{..(f  *  r)  -  cj  ,l;  . 
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Let  W  «  Wj  +  W.j  where  W  is  the  total  defornutioii  energy  of  the  system. 
I'he  defornuition  of  symmetrical  with  respect  to  J  .  Hence  at  tlie 

detornwition  there  is  a  horizontal  tangent  J  .  This  is  expressed  as 
3W 

a  =  ■—  =  0  ,  where  a  represents  the  angle  of  rotation  of  the  tangent  to  J  . 
Likewise,  the  horizontal  displacement  of  J  is  zero,  from  wliicli  we  liave 

Thus  we  obtain  tlie  following  system:- 


|[jc  -  2L  (4  .  h)]  L,  -  ,L;|  ,  {2CL  -  FL^  -  f 

[t  *  '■)  ■  '  i  "-ij  -  " 


which  allows  the  calculation  of  C  and  c  in  terms  of  F  ,  Now  F  is 
fictional,  and  hence  equal  to  zero:- 


=  lahlLJjljL 


We  can  also  determine  6,  -  — 

J  i5F 


Translator's  note  This  equation  cannot  be  read  from  the  copy  supplied 
for  tr.inslat  ion  ;U  Che  bottom  of  page  38. 


We  now  calculate  the  beam  deflections: 


For 


E,i,y';u)  =  -  (x) 


,2  I  2 

EjIiV^fx)  =  ^  (3  -  B  -  l2D)x^  'I  • 
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For  O^J 


and 


(x) 


Ia  ~  ^  ^  (3  -  B)x^  +  (B  -  l)x  . 

ions  of  the  arched  beam  we  put:- 
tiy , 


2A  6  ■  12  ■  "'  6 

In  order  to  study  the  vibrat 


y,(x)  -  6fj(x)  with  =  6fj(x) 


y^Cx)  =  6f  (x) 


^^2 

dt 


=  6f  fx)  . 


The  kinetic  energy  of  a  section  of  the  beam  of  length  dx  is:- 


2  U) 

and  the  total  kinetic  energy  of  the  whole  arched 


beam  will  be ; 


T  = 


I 


L 


.1 


-2. 


i'_,(x)dx  +  2  J  fj(x)dx 
0  0 


(99) 


The  sum  of  the  total  deformation  energy  W  =  W,  .  W,  and  the  total  kinetic 
energy  T  is  constant,  since  the  system  is  not  dissipati^. 


Hence  W  +  T  = 


constant  and  we  shall  calculate  f^Cx)  and  f^Cx)  . 


r 

f  Ivl  =  _ LA.  3  -  B  -  I2D  2  D  3 

1  2  A -  ^  X 


AE.IjL 


f2(x)  = 


4E2I, 

X 

AE  ,  I , 

L.^ 

IT  3 
2L  X 
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wi-  .sli.ill  usi-  till'  lollowi 


On  ord.-r  to  oahuh.to  tl.o  cool  1  i  o  i  ont  s  A,  IJ  .„ul  1) 
nuHUTiial  valiu-s:- 


"  N/nr  ,  1  =  14^ 


10 


'> .  <>  I  V  10 


- 1  2 


I  I'',.  -  10*^'  r,'^  ■ 


.  I.  -  V  10 


-  i 


111  , 


t  1(1111  will  ell  wt‘  ohtain:- 


lliiis  wo  olitain:- 


A  =• 

II  =  =  \ 

l>  “  I.  !')')  V  !()“'’ 


lT(x)(lx  °  4.7')ti— i 


/ 


l^■;(x)^lx  =  4.2H4  X  10‘^i, 


1(1:- 


Sp  iS  ■ 


'i,  7S(1 


•1.284  \  10^1 


-  }nS  ‘ 


‘Si mi  l.irl  v:- 
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By  replacing  6  by  its  value  in  Wj  and  we  liave:- 


and  since  W  +  T  =  constant  we  tind:- 


iS  ^  4  =  0  ,  where  =  —  .  (101) 

b  b 


W  = 


1,5 


0.  )h 


1 ,251,. 


2  2 


'  1  S  / 


i;  1  V 


Hence 


■> 


11,5 

f0.3h  '-''‘‘A 

.*^2'2  ‘4',) 

i 

Sp  j 

/  8  L, 

4.284  X  10  1.  +  4.75b  — 

[  1. 

4  ) 

Calculation  gives 

4^  =  21bb  Hz  . 

2 II 

Consequently,  if  the  arched  beam  is  excited  at  this  frequency  it  will  start  to 
resonate.  Wc  have  checkeil  this  by  experiment  on  .i  model  which  was  dynamically 
similar  to  the  resonator  (Kig  15). 

The  main  resonance  appears  at  about  2000  Hz  as  showit  in  I'ig  3b. 

The  method  used  to  obtain  this  freciueiicy  is  based  on  the  hypothesis  of  a 
small,  static  deformatioit,  which  is  justilied  here.  On  the  other  hand,  the 
hypothesis  of  a  purely-solid  fixing  foi  the  strips  is  much  more  difficult  to 
satisfy  in  practice.  This  would  explain  t  hi'  ilitferences  between  the  theoretical 
and  experimental  values. 

In  he  two  cases  studied,  the  overpressure  due  to  the  resonance  vibration 
of  the  supports  is  in  excess  ot  25  dll,  giving  a  eoelticient  ol  18  which  multiplies 
the  acceleration  of  the  applied  vibr.it  ion.  for  .i  sm.ill  .icce  ler.it  ion,  g  for 
example,  some  lb  g  are  obtained  .it  reson.ince,  which  could  le.iil  to  bre.ik.ige  or  .it 
least  a  high  degree  of  degraii.it  ion  in  both  the  supi'erts  and  the  crystal. 

Wc  have,  therefore,  sought  to  obt.iin  .i  type  ol  support  lor  which  the 
response  shows  no  over-pressure  in  the  r.inge  2tl-2tHH1  11.-. 

1 


Our  I'lioico  lias  talleu  ou  four  supports  of  trapozoidal  profilo,  hold  hy 
ihoriuii-comprossod  rivets  (Kig  37).  The  KNSCMH  p i ozoo  1  oc t  r i o  laboratory  have  ooii- 
stiiutod  several  of  these  resonators  which  liave  withstood  sinusoidal  vibrations 
with  hi^h  accelerations  and  have  enahleil  the  measurements  to  lie  coiiijileted. 

■file  vibration  experiments  were  carried-out  on  an  analogue  model  (Ki^;  38) 
which  gave  a  response  which  was  iierfectly  flat  over  the  range  ol  lre<iuencies 
being  consideri-d  (.Fig  3b).  Above  3000  Hz  t  be  exciter  itself  is  resonating  and 
the  response  from  the  model  is  false.  The  slope  ol  the  response  curve  shows  a 
slight  over-pressure,  equal  to  I.*'*,  al  2000  Hz. 

4.111  Cone lus ions 

The  influence  ot  the  resonator  supports  is  not  negligible  since  they 
introduce  natural  modes  of  vibration  which  can  lead  to  destruction  of  the  quartz 
crystal.  The  resonaiot  is  sensitive  to  1  ow-f  ri'i|uenc.y  vibrations  and  the 
supports  themselves  should  not  increase  this  effect.  The  choice  of  four  supiuirts 
seems  capable  of  supplying  the  reiju i remen t s  of  mechanical  strength  in  a  severe 
env i ronment . 

Chapter  5  Influence  of  'white  noise'  vibration  on  the  resonator 

■file  study  is  similar  to  that  for  sinusoidal  vibrations  but,  with  'white 
noise'  the  disturbance  covers  all  frequencies  within  the  chosen  range.  The 
measurement  [irinciple  is  the  same.  The  disturbances  applied  were  in  the  range 
2t)  to  2Ul)0  Hz,  adjustable  in  80  steps,  each  ol  25  Hz.  The  power  spectral  density 
ot  the  acceleration  w.is  constant  but  the  amplitude  of  the  vibration  varied  as  a 
function  of  frequency.  The  experimental  apparatus  was  analogous  to  that  used 
lor  sinusoidal  vibrations  as  were  also  the  methods  of  measurement  and  analysis. 

The  application  of  white  noise  allows  us  to  define  the  law  for  the  effect 
ot  the  disturbance  as  a  function  of  frequency.  The  'response'  of  the  crystal 
to  this  excitation  is  an  overall  one,  whilst  in  the  case  of  sinusoidal  vibration 
it  was  a  discrete  one.  We  shall  see  that  the  results  confirm  those  obtained  for 
sinusoidal  vibrations. 

5 . 1  experimental  results 

Once  again  we  shall  consider  three  particular  positions  for  the  quartz 
crystal.  Fig  40  shows  the  spectrum  for  the  undisturbed  resonator  when  the  white 
noise  (2  x  lo  *  g^/Hz)  vibration  acts  in  the  same  direction  as  the  natural  vibia- 
t ion  of  the  quartz  crystal  (vertical  position  ox'  )  there  are  considerable 
increases  in  amplitudes  for  the  very  low  frequencies  (Fig  41).  The  power  spectral 


A3 


density  of  the  frequency  variations,  ,  satisfies  a  law  in  f  ,  which 

agrees  well  with  a  law  for  amplitude  decrease  as  a  function  of  frequency,  in 

-2  .... 
f  ,  as  we  showed  for  the  sinusoidal  vibrations.  The  same  results  are  true  for 

the  vertical  position  oz'  . 

In  the  horizontal  position,  the  quartz  crystal  is  less  sensitive  to  vibra¬ 
tions.  On  the  other  hand  it  is  true  to  say  that  an  increase  in  the  vibration 
has  less  effect  on  the  crystal  in  this  case  than  in  the  previous  positions. 

-I  2 

Fig  42  shows  the  frequency  fluctuations  for  a  vibration  density  of  4.50  x  10  g /Hz 
(about  30  g,  effectively).  It  is  obvious  that  the  vibration  density  must  be  more 

than  doubled  in  order  to  obtain  an  effect  comparable  with  that  of  Fig  41.  Again 

-2  .  .  .  . 

there  appears  to  be  an  f  law  for  the  variation  of  the  modulation  amplitude  as 

a  function  of  frequency.  This  accuracy  is  important  because  the  slope  was  less 

significant  in  the  case  of  sinusoidal  vibrations. 

5. II  Long-term  effect  of  'white  noise*  vibration  on  the  natural  frequency  of 
quartz 

We  know  that  the  vibration  causes  either  a  frequency  modulation,  if  the 
disturbance  is  sinusoidal,  or  fluctuations  in  the  natural  frequency  of  quartz  if 
there  are  several  modulating  frequencies.  It  can  be  asked  how  the  initial 
frequency  behaves  when  the  disturbance  is  applied  for  a  long  time,  that  is  to  say 
if  there  is  a  displacement  in  the  initial  resonant  frequency  during  the  vibration. 
For  that,  it  is  sufficient  to  close  the  servo  loop  and  to  use  the  method  of 
measurement  developed  for  the  continuous  accelerations.  In  fact,  no  precise  law 
comes  out  of  it  to  indicate  how  this  initial  frequency  develops  during  the  vibra¬ 
tion  period.  It  seems,  then,  that  the  conclusion  can  be  drawn  that  the  integral 
with  time  of  the  effect  of  sinusoidal  or  random  disturbances  on  the  quartz  crystal 
is  zero.  In  the  same  way,  there  appear  to  be  no  sensible  variations  in  the 
resistance  to  motion  of  the  resonator  during  sinusoidal  or  random  vibrations.  On 
the  other  hand,  if  the  frequency  of  the  sinusoidal  excitation  corresponds  to  a 
vibration  frequency  of  the  quartz  crystal  and  support  assembly,  the  voltages  at  the 
resonator  terminals  show  very  considerable  fluctuations. 

All  the  measurements  carried-out  for  a  white-noise  excitation  of  the  resona¬ 
tor  confirms  the  previous  results. 

Conclusions  from  Part  I 

In  this  first  part  we  have  studied  the  resonator  alone,  placed  in  a  trans¬ 
mission  mounting  so  that  its  parameters  could  be  measured.  For  continuous 
acceleration  measurements  show  that  the  frequency  variation  depends  upon  the 
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intensity  and  direction  of  the  applied  acceleration  with  respect  to  the 
ciyst.il  lographic  axes  of  the  crystal.  In  tlie  case  of  sinusoidal  vibrations, 
analysis  shows  that  their  effect  is  to  produce  a  low-frequency  modulation  of  the 
qu.irtz  natural  frequency.  'White  noise'  vibrations  lead  to  fluctuations  in  the 
quart/,  crystal  frequency  whose  amplitude  spectrum  varies  as  f  ^  .  The  crystal 
position  is  still  preponderant  in  the  case  of  vibrations.  A  model  of  a  crystal 
i-.wiunting  is  proposed  which  can  withstand  very  large  continuous  .accelerations 
vlOt)  g)  or  vibrations  of  the  order  of  40  g  within  the  range  of  frequencies  used. 

t’art  2  -  OSCILLATOR  SUBJECTED  TO  ENVIRONMENTAL  CONDITIONS 

In  this  second  part  we  go  on  to  study  the  characteristics  of  an  oscillator 
■. objected  to  the  same  disturbances  as  was  the  resonator,  previously.  We  note 
tiic  main  criteria  which  characterise  the  instability  of  an  oscillator  and  the 
;.cthods  used  in  this  case. 

In  order  to  carry-out  the  manipulations,  we  built  an  oscillator  of  medium 
[uality  in  which  all  the  electronics  were  made  compact  by  moulding  in  resin. 

This  precaution  enabled  all  the  tests  to  be  carried  out  using  the  same 
I  lectronics  and  several  resonators.  In  particular,  and  wherever  possible,  we 
allied  out  the  environmental  tests  either  on  the  electronics  alone  but  associated 
with  a  resonator  or  on  the  complete  oscillator, 

Comparison  between  the  results  will  enable  us  to  show  what  an  important  role 
is  played  by  the  resonator. 

Chapter  I  Notes  on  the  characteristics  of  an  oscillator 
'.1  Representation  of  a  quasi-sinusoidal  signal*^ 

The  voltage  produced  can  be  represented  by:- 

v(t)  =  [Aq  +  a(t)J  cos  j^w^t  +  (ftj 

wliere  a(t)  and  ((i(t)  are  random  functions,  which  vary  slowly  with  respect  to 
cos  u^t  ,  and  which  represent  the  amplitude  and  phase  variations,  respectively. 

The  instantaneous  angular  velocity  is: 

io(t)  =  —  (Wq!  +  ift)  =  (Oq  +  (f^Ct) 

where  >j'°(t)  represent  the  fluctuations  of  the  angul.ir  velocity  about  . 


The  effects  of  the  amplitude  fluctuations 
stable  oscillator,  in  comparison  with  the  phase 


are  negligible,  lor  an  ultra- 
I  luctuations  and  so  we  Ctin  write 


The  stability  of  a  quasi-sinusoidal  signal  will  be  studied  in  the  spectral 
region  (purity)  and  in  the  time  field  (instability  l(t)). 

•  •  H  Stability  in  the  Sjiectral  region 


This  is  characterised  by  the  radio  irequeiuy. 


phase  and  1 requency  spectra. 


Tlius  we  define  the  power  spectral  density  of  the  f 
of  the  random  process  ,Ht)  by  the  Kourier  transfom  of 
function  R°  (t) 


retiuency  variations 
its  autocorrelation 


V<')  ■  j 


R^(t)  -  t)  ,  witli  llu‘  l>ar  i ti.Iic.il i ng  a  moan  valno  in  limo. 

TI.O  „„y.ical  agoooral  aoaaic,  S,”  (O  ,  will,  „„ly  posUio^o  .alaoa, 

can  be  written:- 


(f) 


/ 


K  U)  cos  2Tii'i  th 


The  stationary  phase  fluctuations  can  also  be  written  as 


S  (f) 

T 


f  „  ,  a  -  i  2  n  f  I  , 

I 


The  phase  and  frequency  spectra  are  related  by 
S°  (f)  =  An'^f'^S^ff)  . 
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It  is  possible  to  introduce  the  reduced  variable 


y(t-)  = 


Hi) 


and  define 


S  (f) 

y 


Mf) 


From  the  experimental  point  of  view,  the  phase  spectrum  S  (f)  is  obtained 

4’ 

by  analysis  of  the  output  voltage  from  an  ideal  phasemeter  whose  input  is  the 
signal  under  study. 

The  spectrum  S  (f)  arises  by  analysis  of  tlie  output  voltage  from  a  perfect 
4> 

frequency  discriminator  to  wtiose  input  is  applied  the  signal  under  study. 


Radio-frequency  spectrum 

The  power  spectral  density  the  signal  vft)  gives  the  power 

distribution  in  the  spectral  region.  This  spectrum  is  obtained  by  direct  analysis 
of  the  signal  from  the  oscillator: 

to 

S  (f)  =  [  R  Ct)e 

V  I  V 


For  ultra-stable  oscillators,  the  phase  and  radio-frequency  spectra  are 
related  by  the  following  equation  when  the  approximation  to  a  small  modulation 
index  is  satisfied: 


S  (f)  -  7 

V  2 


^  V  ^  -  V 


l.lll  Stability  in  the  time  field 

The  instability  in  the  relative  frequency  I(t)  is  characterised  by  the 
standard  deviation  of  the  frequency  variations  (}°(t)  measured  over  a  period  t 


where  <iF°)t,r  is  a  stationary,  random  variable  wliich  represents  the  mean  value 
of  the  oscillator  pulse  rate  over  a  period  of  observation  x  ,  and  hence:- 
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and 


c  2 


and  wc  obtain: - 


o  |<  ji°  )  t ,  T  I 

■"o 


)  ** 

Ot' 

rr  [  S  °(m) 
2-n  J  (f 

fsrn  — 1 

1  du) 

a»T 

1  T  1 

The  importance  of  S ° tm)  shows  up  in  Lliis  equation  since  the  instability 

4> 

Ifi)  can  be  determined  from  S  (uj)  . 

'P 

This  equation  allows  us  to  change  from  the  representation  of  the  stability 
in  the  spectral  region  to  the  representation  in  the  time  field. 

In  practice,  estimations  of  1(t)  can  only  be  made  from  a  large  number  of 
measurements,  independent  of  the  frequency,  over  a  time  period  t  ,  The  character¬ 
istics  in  the  time  field  are  easily  obtained  by  using  time  or  frequency  counters. 
This  possibility  has  led,  in  addition,  to  theoretical  justifications  by  reconinend- 
ing  the  use  of  the  Allan  variance. 

Allan  variance 


If  we  have  a  system  of  N  measurements  of  tlie  relative  frequency,  eacli  of 
duration  t 


h*’ 


'■k  ■  T  / 


ylO)  dt' 


where  k  =  1 ,  . . . .  N 

^.1  =  ^  ^  T  . 

where  T  is  the  period  between  the  start  of  two  successive  measurements,  t  is 
the  measurement  period. 

The  variance  of  this  system  of  measurements  is:- 


c’y(N,T,i) 


>'k 


l.T 
1  Sd; 


which  is  a  random  variable  which  can  be  defined  by  its  mean  <o^(N,T,t))  .  Under 

the  best  case,  (o^(N,T,t))  is  very  slowly  convergent  wlien  N  ‘  .  Very  often 

2 

<0  (N,T,t))  diverges  (when  N  ‘  '■»’)  and  in  practice  it  is  very  difficult  to  allow 
N  to  increase  indefinitely.  Thus  N  has  to  be  limited  and  adjacent  measurements 
have  to  be  imide.  The  Allan  variance  currently  used  is  defined  by:- 


o^(2,i , t) 


denoted 


and  corresponds  to  the  ideal  variance  calculated  from  an  infinite  luunber  of  pairs 
of  measurements. 

The  experimental  estimate  gives:- 


Oy(t) 


1  V  ~ 

M  L, 

k=l 


This  variance  has  the  advantage  that  it  converges  for  the  types  of 
noise  applied,  here,  to  the  oscillators. 

I .V  Noise  in  the  oscillators 

Fig  43  is  a  schematic  diagram  of  the  oscillator.  The  noise  arising  inside 
the  oscillation  loop  is  called  'internal';  that  of  the  output  circuits  is  called 
'external ' . 

The  macroscopic  noise  is  due  to  temperature  variations,  to  shocks  and  to 
variations  in  the  supply  voltage.  Their  effect  is  most  noticeable  in  the  study 
of  the  long-term  stability. 

The  microscopic  noise  is  due  to  thermal  electron  agitation  and  to  noise  in 
the  semi-conductors  themselves. 

Internal  additive  thermal  noise 

This  arises  in  the  oscillating  loop  and  forms  tlie  ultimate  natural  disturb¬ 
ance  suffered  by  the  oscillator.  Numerous  studies  have  led  to  the  following 
equations :- 


50 


I(t) 


where  k  is 
T  is 
P  is 
0  is 
T  is 

The  stability  curve  for  I(t)  has  a  slope  which  varies  at  t  ^  . 

External  additive  thermal  noise 

The  effects  of  this  noise  are  studied  by  superimposing  a  noise  voltage 
e(t)  ,  of  thermal  origin,  onto  the  signal  from  the  loop  (and  supposed  to  be  noise- 
free).  This  noise  voltage  arises  in  the  circuits  external  to  the  loop. 


the  Boltzmann  constant 

the  effective  temperature  of  the  noise  source  (K) 
the  power  delivered  to  tlie  resonator  (W) 
the  selectivity  coefficient 
the  measurement  period  (s). 


There  appears  to  be  a  component 
signal  Aq  ,  which  gives  rise  to  phase 


e^(t)  ,  in  quadrature  with  the  amplitude 
variations;- 


<P  (t) 
e 


e  (t) 
_ 


and  to  frequency  instabili tv:- 


I  (x) 

e 


(OoT 


and  if  u  x  »•  1 

e 


I  (t) 
e 


where  ui 

c 


P 


is  the  filter  cut-off  frequency 

is  the  power  of  the  noise  generated  in  the  output  circuits 

1 

2 


This  noise  shows  a  slope  of  t 


for  the  stability  I(t)  . 
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Scintillation  noise 

The  physical  causes  of  frequency  'flicker'  noise  are  the  slow  and  random 
variations,  of  macroscopic  origin,  of  those  parameters  which  define  the  mean 
frequency.  The  result  is  a  frequency  modulation  and  the  observance  of  a  1/f 
spectrum  at  very  low  frequencies.  This  noise  causes  an  instability  I(t) 
independent  of  t(t  »  1  second),  which  thus  produces  a  limit  in  the  oscillator 
stability. 

l.VI.  Experimental  method  of  measuring  the  instability  I(t) 

Principle  (Fig  4.4) 

A  frequency  comparator  with  error  multiplication  is  used  to  elaborate  a 

frequency  whose  fluctuations  represent  K  times  those  of  the  oscillator  being 

compared.  This  signal,  of  1  MHz  +  KAf  controls  a  'standard  time  interval 

generator 'which  supplies  pulses  separated  by  an  interval  equal  to  10^  periods  of 

the  input  signal.  These  pulses  control  the  starting  and  stopping  of  a  very 

P 

accurate  chronometer,  which  thus  counts  the  duration  corresponding  to  10  periods. 
The  variations  in  this  measurement  allow  us  to  calculate  the  instability  in  the 
time  field  for  a  sampling  period  of  . 

Frequency  comparator  with  error  multiplication 

The  frequency  f^  ,  considered  as  a  reference,  is  multiplied  by  a  factor  of 
10  and  then  mixed  with  the  frequency  fj  multiplied  by  a  factor  of  9.  Each  stage 
elaborates  this  operation.  The  output  signal  F(t)  ,  therefore  contains  the 
frequency  fluctuations  multiplied  by  the  coefficients  of  four  stages. 

The  pass-band  of  the  comparator  input  circuits  is  about  4  kHz.  However, 
between  stages  2  and  3,  and  between  3  and  4,  there  are  two  quartz  filters  with 
narrow  pass-bands  (about  50  Hz) . 

The  consequences  of  this  band  width  can  appear  in  the  calculation  of 

stability:-  B  =  w  -  u 
H  B 
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To  obtain  the  natural  noise  of  the  measurement  chain  it  is  sufficient  to 
apply  the  same  signal  to  the  two  circuits. 

This  system  ot  measurement  appears  to  be  sufficient  to  study  oscillator 
drifts,  but  its  short-term  characteristics  limit  it  to  the  study  of  medium- 
quality  oscillators. 

In  the  following  paragraphs  we  shall  see  that  another  possibility  exists 
for  measuring  the  stability  of  an  oscillator  when  the  s;unpling  time  is  1  second. 

Second  method  o\  miasuring  1 t i )  when  t  =  I  second  -  drift  measurement 

A  second  method  used  ti>  mi'.isure  the  stability  of  an  oscillator  over  I  second, 
as  well  as  to  evaluate  its  dritt,  consists  of  dividing  tlie  reference  and  oscillator 
signals  to  obtain  a  train  ol  111?,  pulses.  This  operation  is  carried  out  using 
fixed  dividers  tN  “  '  10  1  ol  high  quality  and  of  high  accur.icy.  Kach  train  of 

pulses  initiates  the  starting  and  stopping  of  a  very  accurate  chronometer  througli 
the  intermediary  of  a  chronology  control.  The  accuracy  of  this  method  is  limited 
only  by  fluctuations  in  the  reterence  tatomic  clock)  and  the  resolution  of  the 
chroiKvueter  (,'1  ..s). 

The  chronology  control  ensures  that  orders  are  given  to  reset  to  zero,  to 
place  in  the  meuK’ty  and  to  record.  Thus,  every  second,  we  obtain  the  measurement 
ot  the  error,  0  ,  between  the  trains  of  pulses  from  the  reference  and  from  the 
oscillator  under  test  tFig  •451. 

Comparison  chain 

Fig  4b  shows  a  block-diagram  of  tlie  system.  All  the  measurements  are 
recorded  on  nugnetic  tape  to  a  given  format.  The  sampling  time  t^^  determines 
the  rhytlun  by  which  the  measurements  are  obtaiiaed.  Kach  measurement  is  labelled 
with  the  time  tminute,  second,  tenth  of  second)  at  which  it  was  nude.  We  also 
know  the  time  history  of  the  measurements,  and  this  enables  the  oscillator  drift  to 
be  known. 

Recording 

Kach  measurement  'word'  is  formed  from  25  cliaracters  of  which  two  indicate 
the  start  and  end  ot  a  word.  Tlie  word  contains  the  time  as  well  as  the  measurement 
from  the  two  chronometers  used.  Individual  recording  blocks  are  formed  from 

words  separated  by  an  indication  of  tlie  end  of  the  block.  Kach  test,  containing 
a  certain  number  of  blocks,  const i tutes  one  set  of  data  which  is  distinguished  by 
the  appropriate  orders. 
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Depending  upon  Lhe  sampling  time,  ,  the  time  taken  to  record  a  block 

varies  from  20  seconds  to  3  minutes. 

Data  analysis 

The  recordings  are  analysed  witli  a  program  on  an  IRIS  80  computer.  However, 
between  recording  and  data  analysis  there  takes  place  an  operation  which  trans¬ 
forms  the  original  tape  (.7— channel  LBCDIC)  into  one  compatible  with  the  9— channel 
EBCDIC  reading-head  ol  the  computer. 

The  calculations  carried-out  give  both  the  instability  1(t)  ,  calculated 
from  the  Allan  variance,  and  also  the  drift  corresponding  to  a  measurement  period 
of  one  or  several  recorded  blocks. 

Chapter  2  Oscillator  subjected  to  a  continuous  acceleration 

We  use  the  two  methods  of  measurement  of  the  instability  I(t)  which  we 
have  described  earlier.  The  first  gives  Itt)  for  values  of  t  between  1  ms 
and  I  second,  and  the  second  allows  the  oscillator  drift  to  be  found. 

2.1  Measurement  of  instability,  1(.t) 

1  8 

Analysis  of  the  operation  of  the  error  multiplication  comparator  shows 
that  the  instantaneous  frequency  which  it  gives  out  is  equal  to:- 

Flt)  =  f,  +  K°  U) 


where  f^  is  the  nominal  output  frequency 
and  K  is  the  multiplication  coefficient. 

From  this  the  total  phase  of  the  signal  F(t)  can  be  deduced, 


-fiO  =  Oj^^t  +  KvHt) 


The  phase  variation  during  the  time  interval  U.t  +  t)  is  equal  to; 


(f(t  +  t)  -  f)(,t)  =  U)pT  +  Rlf'it  +  t)  -  if(t)] 


and  hence:- 


T  +  — (ifit  +  t)  -  $U)1 

“o 
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Ni>w  t  lu'  moan  valuo  ol  ilu'  I  ri'inn-noy  varialions  liuriii}’,  tlio  l  inu- 
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J.Il  I'.xpo l  imontal  la’sul  Is 

Wo  liavo  I  raoi-il  tlio  varialiiiiis  ol  l(i)  as  a  Iniuilioii  ol  l  lu'  appliod 
ai‘ 1' I  ora  t  i  on  wlu'ii  i  I  aki's  ililloront  valiu's  O' i ‘W)  .  It  i-an  hi'  soon  that  t  ho 
slopos  ol  t  ho  onrvos  liiiiiiiiish  as  i  ilooroasos.  It  would  api'oao,  Ihoii,  that  t  ho 
short-tonii  instaliility  is  dop.railoil  loss  ipiiokly  hy  a  oontiinioiis  aooo  1  ora  t  i  on  than 
is  t  ho  mi'd  intii- 1  o  nil  i  ns  t  ah  i  1  i  t  y  .  This  oont  i  riiis  t  ho  rosiilts  t  roiii  I  ho  rosoiiator 
alono  whioh  p,ivo  tho  doviation  ol  tho  I  roiinonoy  I  roiii  I  ho  iinartz  orystal  iiiulor  tlio 
I'lloot  ol  an  oso i 1  I  a t i on .  This  doviation  oor rosponds ,  horo,  to  an  ' aooo 1  oral od ' 
drill  ol  tho  osoillator  Iroiiiioiioy.  lloiioi'  tho  I  oii>’,-t  oriii  stability  is  dopradod 
iiHiro  qiiiokly  by  tho  aoooli'iat  ion. 


Tiy,  Ab  shows  tho  shapos  ol  tlio  onrvos  ol  Hi)  lor  a  partionlar  position 
ol  tho  quart/,  orystal  and  lor  dilloront  valui's  ol  tho  aooo  I  ora  I  i  on . 


I'ip,  A'l  shows  tho  inllnonoo  ol  tho  orystal  position  with  rospoot  to  l  ho 
aooolorat  ion.  Tho  iiioasnroiiii'iU  t  iiiio  i  »  U)0  ms  is  tho  sanio  Uir  all  tho  onrvos. 
Tho  inoroaso  in  I t i 1  as  a  Innotion  ol  aoooloralion  is  soon.  ITio  variations  ol 
1  (.  i )  do  not  show  tho  saiiio  sliqios  and  this  ditloroiioo  oaii  bo  a  I  t  r  i  but  I'd ,  at  loast 
in  part,  to  tho  iiiothod  ot  mount  inp,  tho  orystal. 
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2. Ill  Oscillator  drift 

The  effect  of  the  continuous  acceleration  is  to  cause  a  displacement  in  the 
'phase-frequency'  curve  for  the  resonator.  An  'accelerated'  increase  of  the 
oscillator  drift  might  thus  be  expected,  liue  to  the  variations  in  the  resonant 
frequency  of  the  quartz  crystal.  Measurements  of  tlie  error  0  ,  made  as  described 
earlier  (Fig  46),  allow  us  to  know  the  variation  of  0  over  a  known  period  as 
well  as  its  sign  with  respect  to  the  initial  error,  0^^  ,  corresponding  to  the 
oscillator  at  rest. 

As  far  as  drift  is  concerned,  we  can  let 

y(t)  =  +  (’,t  +  n(t) 

where  y(t)  is  tlie  instantaneous  frequency  deviation 
Cy  is  the  frequency  at  t  =  t^ 

C|  is  the  slope  of  the  drift 
and  n(t)  is  the  random  noise  in  the  oscillator. 

For  any  one  position  of  tlie  quartz  crystal  and  for  a  given  acceleration, 
observation  of  the  variations  of  0  over  a  period  At  allows  C|(l')  to  be 
determined.  The  figure  gives  the  variations  in  0  (observed  every  second  over 
a  period  of  120  seconds)  as  a  function  of  the  acceleration.  Fig  45  shows  that 
if  the  oscillator  frequency  increases,  0  decreases,  and  vice-versa,  which 
leads  to:- 


C|  (I  ) 

>  0 

if 

AO 

<  0 

c,  (r) 

<  0 

if 

AO 

"  0 

Fig  50  summarises  the  results.  The  increases  in  0  for  positions  (1), 

(3)  and  (5)  of  the  crystal  should  be  noted.  These  correspond  very  well  with  a 
decrease  in  the  resonator  frequency  (Fig  7)  and  hence  of  that  of  the  oscillator. 

Positions  (2),  (4)  and  (6),  for  which  0  is  negative,  do,  in  fact,  indicate 
an  increase  in  the  oscillator  frequency  and  confirm  the  results  found  for  the 
resonator  alone. 

Wc  have  tried  to  apply  the  acceleration  to  the  electronics  alone,  whilst 

maintaining  the  operation  of  the  oscillator.  However,  in  order  to  do  this  the 

resonator  must  be  placed  at  the  centre  of  a  rotating  arm  and  the  electronics  at 

LT  its  extremity.  This  requires  the  use  of  long  lengths  of  cable  which  create 
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variable,  parasitic  capacities  and  upset  the  oscillator  frequency.  It  is  not 
then  possible  to  use  the  comparison  chain  for  the  measurement  of  Ifti  , 

Analysis  of  phase  noise 

Fluctuations  in  the  oscillator  phase  have  been  analysed  by  a  method 
described  in  tlie  next  chapter.  The  analysis  shows  that  there  is  only  a  slight 
increase  in  phase  noise  for  the  low  frequency  components. 

Chapter  3  Oscillator  subjected  to  a  sinusoidal  acceleration 

Ihe  same  vibrations  which  were  applied,  previously,  to  the  resonator  alone 
were  applied  to  the  oscillator.  All  the  measurements  carried-out  will  enable  tlie 
various  results  to  be  compared  with  tliose  for  tlie  resonator.  In  order  to  deter¬ 
mine  the  influence  of  the  external  sinusoidal  vibration  in  the  sjiectral  region 
we  have  carried  out  an  analysis  of  tlie  oscillator  phase  noise. 

1 9 

3.1  Measurement  ol  phase  noise  S  (Fig  51) 

- ^ ^ _ i _ 

The  reference  oscillator  (synthesiser)  was  scrvoed  in  phase  with  the 
oscillator  under  test.  For  that  we  used  a  mixer  which  gave  a  voltage  proportional 
to  the  phase  fluctuations  outside  the  pass-band  of  the  servo  system.  The  time- 
constant  (several  seconds)  reduced  the  servo  system  pass-band  to  very  low 
frequencies.  A  low-pass  filter  (0  to  500  kHz)  eliminated  the  possible  component 
at  the  oscillator  trequency.  A  low  noise  amplifier,  of  gain  1200,  gave  direct 
access  for  the  analysis  of  the  signal  in  the  range  0  to  5000  Hz  with  a  lOHz  filter. 
It  was  also  possible  to  record  the  noise  and  to  analyse  it  by  means  of  a  program¬ 
mable  spectral  analyser.  Both  methods  gave  the  same  results  and  thus  the  phase 
spectrum  S  was  obtained. 

3. I I  Kxperimental  results 

The  oscillator  was  subjected,  in  succession,  to  the  same  sinusoidal  vibra¬ 
tions  as  was  the  resonator.  The  very  low  frequencies  appeared  very  clearly  in 
the  phase  spectrum.  With  respect  to  the  S  spectrum  when  the  oscillation  is 
not  disturbed  (Fig  52),  it  can  be  seen  that  the  80  Hz,  10  g  (Fig  53)  vibration 
not  only  appears,  but  that  it  increases  the  noise  level  by  about  30  dB.  A  380  Hz, 
20  g  vibration  (Fig  54),  shows  a  smaller  effect  on  the  oscillator  phase  noise. 

When  the  vibration  frequency  is  increased,  the  disturbance  peaks  are  again  present 
in  the  oscillator  phase  spectrum.  Figs  55,  5b  and  57  show  the  effect  of 
frequencies  equal  to  1750  Hz,  2700  Hz  and  4080  Hz,  respectively,  all  for  the  same 
20  g  excitation.  For  a  frequency  of  around  4000  Hz,  the  natural  resonance  of  the 
exciter  intervenes  to  generate  other  f retpieneies,  which  explains  the  other  peaks 
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close  to  the  excitation  Irequency,  The  decrease  in  amplitude  of  these  different 
peaks  should  be  noted. 

In  order  to  obtain  a  direct  and  continuous  measurement  of  tlie  effect  of  the 
external  vibration  f requeitcy  we  used  the  system  shown  in  Kig  58. 

The  vibration  generator  was  controlled  in  frequency  by  tlie  spectrum  analyser 
itself  which  displaced  its  analysis  filter  at  a  given  velocity.  There  was  also 
synchronisation  between  the  external  vibration  frequency  and  the  central  analysis 
frequency.  Hence  we  were  able  to  obtain  the  oscillator  phase  noise  and  its 
response  to  a  sinusoidal  vibration  in  the  range  20-2000  Hz  (Kig  59). 

3 

The  decrease  in  noise  obeys  a  1/f  law.  The  position  of  the  crystal 
the  disturbance  is  also  critical. 

It  is  interesting  to  note  that  tl»e  vibration  frequencies  appear,  here,  very 
clearly  superimposed  on  tlie  natural  noise  of  the  oscillator.  In  the  study  of 
the  resonator,  tlie  analysis  did  not  permit  frequencies  above  500  Hz  to  be  detected. 
In  the  present  case  we  see  that  the  crystal  is  effectively  sensitive  to  vibra¬ 
tions  greater  than  500  Hz  and  their  detection  is  facilitated  by  the  operation  of 
the  oscillator. 

It  can  be  stated  that  the  oscillator,  disturbed  in  this  way,  sees  its 
frequency  modulated  at  the  frequency  of  the  external  vibration. 

Chapter  A  Oscillator  subjected  to  a  pseudo-random  'white  noise'  vibration 

The  external  vibration  applied  to  the  oscillator  is  contained  within  the 
range  between  20  and  2000  Hz.  The  oscillator  used  was  the  same  on  as  was  used 
previously  and  the  resonators  used  were  those  already  studied  in  I’art  1. 

Analysis  of  applied  noise 

The  vibrations  were  maintained  at  a  constant  power  spectral  density  within 
the  range  of  frequencies  chosen.  Spectral  analysis  of  the  applied  noise  showed 
a  climbing  of  the  spectrum  beyond  2000  Hz  and  some  peaks  about  3500  Hz.  These 
latter  frequencies  correspond  to  the  resonance  vibrations  of  the  exciting  element. 

4 . I  Instability  I(t) 

Kig  60  shows  a  comparison  of  the  results  obtained  when  the  oscillator  is 

3  2 

subjected,  successively,  to  3  x  I0  g  /Hz  at  2.5  g  ,  5  g  ...  and  10  g  ....  for 

e  1  I  e  t  1  e  1 1 

the  three  usual  positions.  The  arrangement  which  corresponds  to  the  horizontal 
position  of  the  quartz  crystal  has  an  instability  curve  whose  slope  is  always  in 
T  '  and  values  greater  than  those  for  which  the  oscillator  is  not  disturbed. 

The  increase  in  the  disturbance  seems  to  enhance  the  stability  degradation. 
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The  two  other  vertical  positions  of  the  quartz  crystal  show  a  very  high 

range  of  instability  values  whose  mean  slope  is  sensibly  modified  to  the  form 

-0.6 

t 

This  modification  of  the  slope  is  no  doubt  attributable  to  the  increase  in 
the  internal  natural  noise  of  the  resonator. 

On  the  other  hand  it  would  appear  to  be  difficult  to  relate  each  instability 
value  to  a  level  of  applied  disturbance. 

We  have  measured  the  instability  of  the  same  oscillator  placed  under  the 
same  vibration  conditions  but  with  the  resonator  screened  from  the  disturbance. 

The  measurements  gave  the  instability  curves  shown  in  Fig  61,  for  which  we  see 
a  slight  increase,  but,  however,  less  than  the  smallest  instability  recorded  when 
the  resonator  was  also  disturbed. 


This  instability  appears,  moreover,  to  be  independent  of  the  vibration 
applied  to  the  oscillating  electronics  alone. 

4. II  Oscillator  drift 


Using  the  second  method  of  measuring  I(t)  when  t  =  1  second,  we  calculated 
the  variation  of  the  initial  displacement  as  a  function  of  time  passed  and 
applied  acceleration.  The  results  showed  a  drift  which  was  always  negative  (the 
oscillator  frequency  increased)  and  the  highest  values  were  obtained  for  the  two 
vertical  positions  of  the  quartz  crystal. 


4. Ill  Phase  spectrum  of  an  oscillator  disturbed  by  'white  noise' 

We  used  the  same  experimental  apparatus  as  has  already  been  described  for 
the  sinusoidal  vibrations.  Fig  62  shows  the  phase  spectrum  S  for  the  undis- 

't’ 

turbed  oscillator.  When  the  quartz  crystal  was  in  the  vertical  position  the 

spectrum  S  was  very  clearly  disturbed  (Fig  63)  within  the  frequency  range 

0  to  2000  Hz.  The  effect  was  slightly  less  for  the  horizontal  position  of  the 

resonator.  The  phase  noise  returned  to  its  normal  level  above  2000  Hz,  and  this 

phase  noise  was  recorded  for  analysis  by  means  of  a  program.  The  slope  of  the 

3 

notse  appeared  very  strongly  (Figs  64  and  65)  to  be  of  the  form  1/f  and  we 
have  identified  it  with  a  frequency  flicker  noise.  The  curves  have  been  super¬ 
imposed  in  order  to  facilitate  comparison  between  the  spectra  when  the  oscillator 
was  undisturbed,  (curves  1  of  Figs  64  and  5  of  Fig  65),  then  disturbed  in  the 
horizontal  position  (5  g  curve  2,  10  g  curve  3)  and,  finally,  in  the  vertical 
position  (5  g  curve  6,  10  g  curve  7). 
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Oscillator  noise  was  measured  when  the  electronics  alone  were  subjected  to 
vibration  without  the  resonator.  This  analysis  (Fig  66)  shows  that  the  vibration 
has  a  slight  effect  on  the  electronics  by  introducing  a  noise  whose  slope  is 
again  of  the  form  1/f^  .  However,  the  noise  level  is  mucli  less  tlian  when  the 
quartz  crystal  was  affected  by  the  external  vibration  along  with  the  oscillator. 

In  particular,  for  the  oscillator  alone,  the  phase  noise  returned  to  the 
normal  level  above  1000  Hz  (with  vibration  applied  within  the  range  20-2000  Hz), 
whilst  in  the  case  where  the  oscillator  was  complete,  the  applied  white  noise 
acted,  effectively,  over  the  whole  range,  that  is  to  say  up  to  2000  Hz. 

Conclusions  from  Fart  2 

A  quartz  oscillator  disturbed  by  a  continuous  acceleration  shows  a  long¬ 
term  frequency  variation  which  is  a  function  of  tlio  intensity  and  direction  of 
the  applied  acceleration.  When  subjected  to  a  sinusoidal  vibration,  its  frequency 
is  modulated  by  the  vibration  frequency  found  in  its  phase  spectrum. 

The  effect  of  'white  noise'  is  to  produce  a  pliase  noise  in  the  accelerator 

-3 

of  form  f  ,  that  is  to  say  a  frequency  noise  flicker.  Measurements  of  the 
instability,  1(t)  ,  show  the  variations  as  a  function  of  the  applied  acceleration 
and  of  the  position  of  the  oscillator  with  respect  to  the  direction  of  the 
disturbance . 

It  is  interesting  to  note  the  presence  of  frequency  flicker  noise  when  the 
oscillator  is  excited  with  'wl\ite  noise'.  It  can  be  said,  then,  that  the  random 
vibrations  are  one  of  the  origins  of  tliis  noise,  and  are  very  important  in  the 
study  of  oscillator  stability. 

General  conclusions 

In  the  first  part  of  this  study  we  were  interested  in  quartz  resonators 
alone,  subjected  to  steady-state, sinusoidal  or  random  accelerations.  From  both 
the  theoretical  aspects  and  the  experimental  results  we  have  been  able  to  describe 
the  behaviour  of  the  crystal  in  each  case.  Wlien  it  undergoes  a  steady  accelera¬ 
tion,  six  positions  determine  the  variation  in  the  crystal  resonant  frequency 
perfectly,  and  show  it  to  be  linear,  to  a  first  approximation. 

Sinusoidal  or  random  vibrations  modulate  the  natural  frequency  of  the 

resonator  at  very  low  frequencies,  or  create  fluctuations  in  the  natural 

-4 

frequency  whose  spectral  density  varies  as  f 

The  different  results  arising  from  measurements  on  the  oscillator,  subject 

LT  to  the  same  disturbances,  showed  the  importance  of  tlie  contribution  from  the 
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resonator  in  the  oscillator  frequency  fluctuations.  It  would  appear,  then,  that 
the  frequency  fluctuations  of  a  quartz  oscillator,  when  subjected  to  severe 
environmental  conditions,  are  essentially  due  to  the  effects  of  the  disturbances 
on  the  quartz  resonator. 

In  particular,  random  vibrations  gave  a  phase  spectrum  in  f  ^  ,  that  is 
to  say  a  frequency  noise  flicker  whose  origin  is  shown  in  this  paper. 

These  results  should  have  a  beneficial  result  in  determining  the  best 
conditions  for  the  production  of  quartz  oscillators  designed  for  aerospace 
applications . 

The  next  studies  should  be  orientated  towards  seeking  new  ways  of  cutting 
the  quartz  which  might  decrease  the  effects  of  the  environmental  conditions. 
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Apparatus  for  measuring  the  frequency  variation  of 
quartz  resonator 


Figs  2-5 


Fig  2  Layout  for  measuring  the 
quartz  parameters 


Fig  3  Servo- system 


Fig  4  Phasemeter  output 
voltage 


Fig  5  Quartz  crystal  cut 
along  AT 


Variations  in  resonant  and  anti-resonant  frequencies 
for  acceleration  directions  (1)  and  (2)  (c/  Fig  6) 
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Fig  9  Variations  in  resonant  and  anti -resonant  frequencies 
for  acceleration  directions  (3)  and  (4)  (c'f  Fig  6) 
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Fig  10  Variations  in  resonant  and  anti-resonant  frequencies 
for  acceleration  directions  (5)  and  (6)  {of  Fig  6) 
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f'igs  15&16 


Displacement  of  the  phase-frequency 
curve  ^or  a  quartz  crystal  under  the 
effect  of  an  acceleration 


'Hysteresis'  effect  in  the  frequency 
variation 
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Fig  25  Quartz  crystal  in  the  horizontal  position  30  Hz  -  4  c 


Figs  30&31 


Fig  30  Mounting  analogous  to  that 
of  the  resonator 


Key; 

amplificateur 

de  charge  =  load  amplifier 

enregistreur 

graphique  =  graph-plotter 

Fig  31  Apparatus  for  measuring 
the  vibrations 
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Fig  32  Response  curve  for  the  support-quartz  crystal  assembly 


Fig  33  Symmetrical  vibrations 
of  an  arched  structure 


Fig  35  Mounting  analogous  to 
that  of  the  resonator 


Fig  34  Analysis  of  the  forces  applied 
to  the  arched  structure 


Fig  40 
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Fig  40  S  -  for  an  undisturbed  quartz  crystal 
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Fig  41 
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Fig  41  S  ,  for  a  quartz  crystal  in  the  vertical  position 

ATq 

^  -12 
subjected  to  white  noise  of  2.10  g  /Hz  within  the 

range  20-2000  Hz 
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Fig  42 


S  ^  for  a  quartz  crystal  in  the  horizontal  position 

Q 


subjected  to  white  noise  of  4.5  x  10'^  g^/Hz  within 


the  range  20-2000  Hz 
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Figs  43&44 
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Fig  44  Apparatus  for  measuring  I(t) 
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Fig  46  Block  diagram  of  system  for  drift 
measurement 
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Figs  50451 


Fig  50  A9  as  a  function  of  r 
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Fig  51  Equipment  for  measuring  S 
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Key; 
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Fig  58  Equipment  for  meaf^i’^i^The  influence  of 
sinusoidal  disturbance  on  the  oscillator 
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Fi9  66  for  an  oscillator  of  which  only  the  electronics  are 
disturbed  by  white  noise 


